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ABSTRACT 


An instrument, the Caltech High Energy Isotope 

Spectrometer Telescope (HEIST), has been developed to 
measure isotopic abundances of cosmic ray nuclei in the 
charge range the energy range between 30 and 

800 MeV/nuc by employing an energy loss -- residual energy 
technique. Measurements of particle trajectories and 
energy losses are made using a multiwire proportional 
counter hodoscope and a stack of Csl(TI) crystal 

scintillators, respectively. A detailed analysis has been 
made of the mass resolution capabilities of this 
I nstrument . 

Landau fluctuations set a fundamental limit on the 
attainable mass resolution, which for this instrument ' 
ranges between 'v.Qt AMU for Z 3 sbd '^•2 AMU for Z-s^8b- 
Contributions to the mass resolution due to uncertainties 
in measuring the path-length and energy losses of the 
detected particles are shown to degrade the overall mass 
resolution to between *^.1 AMU ( Z 3 ) and ^.3 AMU 
(Z-N>8b) • 

i 

A 

A formalism, based on the leaky box model of cosmic 
ray propagation, is developed for obtaining isotopic 
abundance ratios at the cosmic ray sources from abundances 
measured in local interstellar space for elements having 
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three or more stable isotopes, one of which is believed to 
be absent at the cosmic ray sources. This purely 
secondary isotope is used as a tracer of secondary 
production during propagation. This technique is 
illustrated for the isotopes of the elements 0, Ne, S, Ar 
and Ca. 

The uncertainties in the derived source ratios due to 
errors in fragmentation and total inelastic cross 
sections, in observed spectral shapes, and in measured 
abundances are evaluated. It is shown that the dominant 
sources of uncertainty are uncorrelated errors in the 
fragmentation cross sections and statistical uncertainties 
in measuring local interstellar abundances. 

These results are applied to estimate the extent to 
which uncertainties must be reduced in order to 
distinguish between cosmic ray production in a solar-like 
environment and in various environments with greater 


neutron enrichments. 
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i. INTRODUCTION 

Measurements of the nuclear composition (both 
elemental and isotopic) of material on the earth, in 
meteorites, in the sun and stars, and in the interstellar 
medium have contributed greatly to our understanding of 
the astrophy s I ca I conditions under which this material was 

I 

synthesized. In fact, such measurements have been 
essential to the development of the theory of 


nucleosynthesis. 

The 

material observed i 

n these 

s i tes 

generally exists 

I n 

therma I equ i I 

I br I urn 

with 

its 

Surroundings and 

th I s 

has probably 

been 

the case 

ever 

since it was synthesized 

• 




The nuclei in 

the cosmic radiation 

form 

a samp I 

e of 


matter which can be observed at kinetic energies greatly 
In excess of thermal energies (cosmic rays have been 
observed at energies up to 10^^ eV (Brownlee et a!., 

At present it is not clear whether this material 
IS formed in sites different from those which produce the 
matter observed in the solar system, or whether the same 
sites produce both populations of matter with the cosmic 
rays simply being that small fraction which subsequently 
gets accelerated to near-re I at i v i st i c energies. Not only 
IS it possible that the cosmic rays will provide 
information concerning astrophys i ca I objects other than 
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those which produced the solar system material, but they 
also may contain information concerning the time evolution 
of such sources. Measurements of radioactive nuclides in 
the cosmic rays (Hagen et al., Garcia-Munoz et al., 

X^llb; Webber et al., I'll!) lead to estimates of 
~i0 years as the age of the cosmic rays. This is 
much younger than the age of 'solar system material, which 
has been estimated to be years. 

Improvements in the instrumentation used for 
measureme,nt 5 of the nuclear component of the cosmic rays 
have led to measurements of X) the abundances of 
individual elements in the cosmic rays in the charge range 
isotopic abundances of the elements H and He, 3) 
elemental abundances in the charge range Z > H) mean 

masses of elements with and 5) isotopic 

abundances of elements with 3 ~ Z 8 • Reviews of these 
improving measurements of the elemental and isotopic 
composition of the cosmic rays have been presented by 
Gar c I a-Munoz (XS13)» Stone (XS13^» Lund (X815)» Meyer 
(iS15) and Waddington (XR11^* In addition, a number of 
instruments which should be capable of isotopic abundance 
measurements for all elements up through Ni have recently 
begun making measurements or will soon do so. It is 
believed that such measurements of the isotopic 

composition of the cosmic rays wilt further our 
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understanding of the astrophys i ca 1 processes Involved in 
the synthesis and acceleration of this material. 

In section 2 of this thesis we discuss the 
development of one such instrument and present a detailed 
analysis of the influence of a variety of instrumental 
parameters on the mass resolution attainable with this 
instrument. In addition to evaluating the expected 
performance of our instrument, this analysis serves to 
indicate those aspects of the measurement in which 
improvements in the instrumentation can s i gn i f i c.ant | y 
improve the capabilities of this instrument or of other 
instruments of similar design. 

Techniques which have been used for cosmic ray mass 
determination involve measuring two parameters such as 
total energy, range in a selected absorbing material, 
velocity, magnetic rigidity, specific ionization rate or 
intensity of Cerenkov light emission. Each of these 
quantities can be expressed as a product of a function of 
velocity and a function of the particle>'s charge and/or 
massi The fact that only a limited number of discrete 
values of charge and mass are possible allows one in some 
cases to determine both of these quantities from 
measurements of only two to the parameters listed above. 
Techniques for measuring charge and mass using 



H 


comb I n at I'ons of measurements of total energy, specific 
ionization and Cerenkov emission have been reviewed by 
Stone ( iS1H) • 

The Caltech High Energy Isotope Spectrometer 
Telescope (HEIST) is designed to measure the mass of 
cosmic ray particles in the charge range 3 ^ Z 28 with 
energies between 30 and 800 MeV/nuc (the energy interval 
varying with particle charge within these limits). An 
energy-loss - residual energy technique is employed. This 
technique 'involves measuring the amount of energy that the 
particle loses in passing through an absorber of known 
thickness and measuring the residual energy with which the 
particle emerges from this absorber. This technique has 
been employed by a number of investigators for measuring 
both elemental and isotopic abundances. References to 
publications describing these measurements can be found in 
the reviews listed above. 

Several innovations have been included in the HEIST 
instrument. Csl(Ti) crystal scintillators are used to 
make the energy loss measurements. The use of crystal 
scintillators rather than the commonly used plastic 
scintillators should significantly reduce the degree of 
saturation of the scintillation efficiency when' the 
scintillator is exposed to the heavily ionizing particles 
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of interest here. The effects of scintillator saturation 
on mass resolution are discussed by Webber and Kish 
Also, a stack of eight scintillators of 
progressively increasing thickness have been employed. By 
summing the energy losses in all detectors prior to th,e 
one in which a particle stops we can insure that the 
absorber thickness is a large fraction of the particle's 
range. This technique significantly improves the mass 
resolution over that obtainable using a thin dE/dx 
detector followed by a thick total energy detector. A 
multiwire proportional counter (MWPC) hodoscope is used to 

make precise measurements ( (T" mm) of the particle's 

IvvPt 

trajectory. These measurements make it possible to make 
corrections for the particle’s angle of incidence and for 
variations of thickness and light collection efficiency 
with position In the scintillators. Also a technique 
whereby signals from two different dynodes of each 
photomultiplier tube are analyzed using a logarithmic 
pulse height analyzer allows us to achieve the large 
dynamic range required in order to measure signals 
throughout the charge range of interest while still 
maintaining an adequate pulse height resolution throughout 
that range. 



b 


Instruments such as the one which we describe should 
be capable, even with the limited statistics normally 
obtained in cosmic ray experiments, of determining whether 
local isotopic abundance distributions originate in a 
source which is drastically different from that which 
produced the solar system material. However, it is 
possible that differences between the cosmic ray sources 
and a solar~like source may be reasonably subtle, arising, 
for example, from differences in the state of evolution of 
the seed materia! which provides the fuel for a common 
nucleosynthesis process. If it is necessary to 
distinguish source differences which are not dramatic, 
then it will be necessary to significantly improve the 
statistical accuracy of the cosmic ray measurements and 
also to obtain more precise values of the nuclear cross 
sections required for calculating the secondary 
contributions to the observed fluxes. 

In section 3 o'f this thesis we present a formalism 
which can be used in certain cases for obtaining isotopic 
source abundance ratios and their uncertainties. This 
treatment is based on the 'Meaky-box model" of cosmic ray 
propagation (Gloeckler and Jokipii, i%S» Meneguzzi et 
3l* ini)* We use an isotope which is believed to be 
absent at the cosmic ray sources as a tracer of the 
production of secondary cosmic rays during propagation. 
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Measurements of the abundances of this isotope and of two 
other isotopes of the same element in local interstellar 
space can then be used to derive the relative abundances 
of the latter two isotopes at the cosmic ray sources. 
This type of analysis makes it possible to interpret 
measurements of the isotopic make-up of a single element 
in a way which should not be affected by charge-dependent 
processes which may be involved in the acceleration of the 
cosmic rays (Casse et al., iRlSa)* 

We consider in detail the interpretation of the 
isotopic composition of the elements 0, Ne, S, Ar and Ca. 
These elements should provide significant astrophys i ca I 
information. In particular, they should make possible'the 
determination of the neutron enrichiTient of the cosmic ray 
sources and the state of evolution of the seed material 
which has been processed by these sources. 

For the case of a solar-like cosmic ray source, we 
evaluate the contributions of uncertainties in a variety 
of parameters-- 1 oca I isotopic abundances, composition of 
parent species, nuclear cross sections and spectral 
shapes--to the uncertainty In the calculated source 
abundance ratios.’ In this way we identify dominant 
sources of uncertainty. We find that uncertainties in the 
I nterpretat I on of isotope measurements obtained with 
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instruments now being developed and flown should be 
dominated by two sources of error: 1) by statistical 
errors in the measured isotopic abundances and 2) by 
uncertainties in the fragmentation cross sections 
necessary for evaluating the secondary contributions to 
the observed fluxes. 

We then narrow our consideration to these two sources 
of error and evaluate, for a wider range of observed flux 
ratios, the relative uncertainties which they will produce 
in the derived source abundance ratios. These results are 
obtained using typical estimates of the counting 
statistics and cross section errors which can presently be 
achieved. However, the uncertainties which we obtain can 
easily be scaled for other conditions. The details and 
limitations of such scaling are also examined. 

We have applied our uncertainty estimates to the 
problem of determining the extent to which statistical and 
cross section errors must be reduced in order to 
distinguish between cosmic ray production in a solar~like 
cosmic ray source and production in various non-solar 
environments. In particular, we derive the number of 
cosmic ray events of the elements S and Ca and the level 
of uncertainty in the fragmentation cross sections which 
must be achieved in order to distinguish various levels of 
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neutron enrichment over a solar-1 ike environment- This 
treatment is based on a model for the production of the 
elements S through Ca in an explosive oxygen burning 
process, as presented by Woosley et al. (iS13)» We find 
that present uncertainty levels should permit one to 
distinguish a source environment enriched in neutrons by a 
factor of ^3 over the solar composition by using 
measurements of the isotopic composition of cosmic ray 
sulfur. In the case of the calcium isotopes, it appears 
that reduction of the important uncertainties will be 
required even to distinguish a solar-like source from a 
source having as much as a factor of 5 enhancement of 
neutron excess over the solar value. 

Finally, we make estimates of the magnitude of the 
corrections which must be made for solar modulation 
effects in order to apply this formalism to measurements 
made in the vicinity of the earth. 
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2. INSTRUMENTATION 

2-i Measurement Technique 

The HEIST instrument is designed to determine the 
mass of cosmic ray nuclei impinging upon it by means of an 
energy-loss - residual energy measurement. The basic 
detector geometry is illustrated in figure 1 The energy 
loss in the first detector is referred to as AE and the 
energy deposited in the second detector is referred to as 
E'. The measurement relies upon the fact that while the 
specific ionization of a charged particle depends only on 
its charge and its velocity (or equivalently its energy 
per nucleon), particles of mass numbers A and A+i that 
have the same energy per nucleon wi I I have total kinetic 
energies with a relative difference of 1/A. This 
difference is less than 2? "for the isotopes of nickel, the 
heaviest nuclei which we are trying to identify. 

As a charged heavy particle passes through matter it 
loses energy predominantly by means of collisions with 
atomic electrons in the medium. In the energy range of 
interest here, rwr|-|00(] MeV/nuc, the specific ionization 
IS given by (Barkas and Berger, 1%H» Janni, ISbb) 
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F I GUR E [ 


Schematic illustration of geometry for 
AE-E' measurement of particle mass. 
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In this expression, 

Z is the particle's charge in units of the 

proton charge, i.b02xi0 Coulombs, 

M is the particle’s mass in units of the 

proton mass, mp='138*28 MeV/c^, 

^ IS the density (g/cm^) of the absorber, 

2 IS the atomic number of the absorber, 

IS the atomic weight of the absorber, 

X IS the ionization potential of the 
absorber , 

m IS the electron mass, Q.5il MeV/c®’, 



p IS the velocity of the c(\arged particle in 

units of the velocity of light, 
10 

c=3-00xl0 cm/sec. 



^'LC, is a correction for the fact that the 
Z f 

electrons of the medium 'are bound in atomic 
she I Is, 

IS a polarization correction which accounts 
for the fact that the presence of the heavy 
ion distorts the shape of the electron 
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cloud of near-by atoms in such a way as to 
partially shield these electrons. 

The polarization and shell corrections are small for the 
range of energies of interest and will be neglected below. 
In addition, it should be noted that if the square root in 
the logarithm is neglected (it varies by /«0.i? from unity 
for protons at | GeV), then equation 1 can be written in 
the form: 



Here the function s(E/M) is the specific ionization (or 
mean rate of energy loss per unit of matter traversed) of 
a proton in the material of interest, and is expressed in 
MeV/(g/cm^). The mean range of a particle is obtained 
from the i ntegra I ; 



Here Rp(E/M)sRji (E/M) is the range of a proton of energy 
E/M. Numerous tabulations of R^, vs. E/M are available 
(Barkas and Berger, 1%H; Janni, i%b; Northcliffe and 
Sch i 1 1 I ng, ) . 
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The above discussion assumes that the charge, Z, of 
the ion does not change as the particle slows down. 
However when the ion slows to a point where its velocity 
is comparable to the velocity of atomic electrons in their 
orbits, it becomes possible for the ion to pick up atomic 
electrons' thereby altering Z. For Z~28 this occurs for 
energies in the range E/M^|0 MeV/nuc (Northcliffe and 
Schilling, iSlO)- This effect, although significant 
enough to be considered in the present measurement, will 
for the sake of simplicity be neglected in the following 
discussion. 


The function R, (E/M) can be approximated with 
reasonable accuracy by a power taw in E/M over decades 
variation of E/M. We can write 




where k=.0iH2 nnm of Cs| and a = l.bbH« Figure shows the 
exact form of Rp{E/M) as well as this power law 
approximation. Figure shows the fractional difference 
of these two quantities. By noting that the particle’s 
range and energy after passing through a layer of 

thickness L are R^j^(E/M)-L and E* respectively, we can 


write: 
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FIGU RE 2 

Proton range tn Csl vs. particle 
energy . 

a) Solid curve is exact relation as 

given by Janni dashed curve is 

power law approximation. 

b) Fractional difference between exact 
relation and power law approximation. 
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By eliminating R^j^(E/M) between 2 3 we obtain an 

implicit relation for M in terms of Z and the measurable 
quantities L, and E': ^ 

0=f^(M.AE,E;L)=R,(^>Rp(|V 

In the power law approximation we can explicitly solve for 
M: 



in what fol lows we shal I assume that the particle's 
charge, Z, Is known. This is reasonable since, due to the 
discrete nature of the nuclear charge and mass, both of 
these quantities can be determined from measurements of ^E 
and E*. Experiments employing the AE-E' technique have 
already succeeded in resolving individual elements with 
<see, for example, the review by Garc i a-Munoz , 

iS13). 
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The uncertainties in the masses which we will derive 
can be expressed as 



In this equation, , is the uncertainty in the energy 
loss in the AE detector due to statistical fluctuations in 
the ionization energy loss process, and 

are the uncertainties in our measurements of AE, EJ 
and L respectively. Equation 5 applies only if the 
various contributions to the mass uncertainty are 
statistically independent. This condition is not strictly 
correct for our measurement, but the error made in using 
it should not be significant. Equation can be used to 
obtain the required partial derivatives of M with respect 
to the measured quantities. Table [ summarizes these 
derivatives. The measurement uncertainties will be 
evaluated in the course of the discussion of the 
I nstrument . 

This section will discuss our mass resolution 
requirements and the limits on the mass resolution 
attainable using the AE-E’ technique. It is possible to 
plot lines of constant mass, M, on a AE vs. E’ graph. 
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This formal procedure can be carried out even for 
non-integral values of M. By dividing the plane in 
this way the events corresponding to a particular element 
can be binned to form a mass histogram. In our discussion 
of our ability to obtain isotope ratios from such a 
histogram we will treat the simple case in which only two 
isotopes have non-neg I i g i b I e abundances. If we know that 
the mass peaks have Gaussian shape and if we accurately 
know the mass resolution, , of the instrument then we 
can use the method of least squares to obtain the 
abundance ratio of these two species. In addition, the 
uncertainty in this ratio can be obtained from the 
curvature of the at its minimum (Bevington, i%R; 
Mathews and Walker, 1S1Q)- The details of this procedure 
are discussed in appendix A. 

Figure 3 shows as a function of the abundance ratio, 
r, the number of events required to obtain relative 
uncertainties of iO^, 20/5, 5055 and 100^ in our estimate of 
this ratio. The number of counts is weakly dependent on 
the mass resolution as well. The solid curves in the 
figure are calculated for -0>3 units and the dashed 
curves are calculated for =0-25 units (where the 
separation between the centers of the two Gaussian 
distributions is defined as one unit). In addition, a 
pair of dotted curves are included in the figure to 



FIGI^ 3 


Number of events required to determine 
abundance ratios with uncertainties of 
[{]%, 50^ and iOO^- Curves are 
shown for mass resolutions of 0-35 and 
0-30 times the mean separation between, 
the distributions. Dotted curves 
indicate combinations of N and r which 
yield i and J.0 events of the less 
abundant species. 
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indicate those combinations of r and N which yield i and 
iO counts in the peak having the smaller abundance. It is 
not expected that this type of calculation of OJ! will be 
reliable when there are ^10 counts in the less abundant 
peak for reasons mentioned in appendix A. 

It should be noted that although distributions of 
known shape can be separated by fitting techniques even 
when they have a large degree of overlap (as indicated by 
the weak dependence of on in figure 3)> 't is 
important that we be able to achieve a reasonable degree 
of separation. This is because this measurement is 
subject to various systematic errors which will tend to 
produce non-Gaussian distributions whose precise shape is 
difficult to determine. As an estimate of the required 
mass resolution we have calculated the mass resolution for 
which an inflection point is obtained between two Gaussian 
distributions whose means are separated by one unit, as a 
function of the ratio of the number of counts in the two 
distributions (Stone, 1S13; Hagen, iSlb)* The details of 
this calculation are outlined in appendix A and the 

t 

results are summarized in figure 4- 

If the iron produced in the cosmic ray sources 
consists of a single isotope. A, one expects a ratio 
N ( A-| ) /N ( A ) ;^ 0 • i a'l' earth due to spallation reactions in 
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the interstel lar medium. Thus we need to achieve a mass 

resolution of £0*30 AMU. There are, however, instances 

where the mass peaks which we are attempting to resolve 

are separated by two atomic mass units. For a given 

element nuclear pairing effects (Evans, 1S55; Preston, 

enhance the stabi I i ty of nuclei with even numbers of 

neutrons (N). Hence, there are cases where an odd-N 

isotope lies between two isotopes of interest, and this 

isotope can be neglected either because it is unstable 
39 

(e.g., Ar) or because its production is suppressed in 

the astrophy s i ca 1 environment of interest and it is not 

produced in significant amounts by spallation reactions 

91 

during propagation (e.g., Ca). In such cases we can 
resolve the two even-N nuclei with a mass resolution as 
large as AMU. 

The fundamental limit on the mass resolution 

attainable by the AE-E* technique is due to the 

statistical nature of the Ionization energy loss process. 
This results in Landau f 1 u ct u at i ons--a distribution of 
possible AE values for a monoener get i c beam of identical 
particles incident on the detector. The energy losses of 
interest here are su f f i c'l ent I y large in comparison with 
the energy that the particle can transfer in a single 
collision so that the Landau distribution is accurately 
approximated by its Gaussian limit (Seltzer and Berger, 



i%4; Ross!, iS52), 


gr’- = (a39fcHe/f^ { 


^ECU'KolaiO 


where 

L, = thickness of AE detector (cnn), 

j 

^ (Lorentz factor), 

^ = density of the detector material (g/cm'^), 
^ = mean atomic number of the detector 

matei" i a I , 

^ = mean atomic weight of the detector 

mater i a 1 . 

This expression applies only for L sufficiently small so 
that the particle’s specific ionization does not change 
significantly as it traverses the AE detector. 


The j4E detectors of interest here are thick enough so 
that a correction must be made for the deceleration of the 
particle. We introduce a "deceleration factor", D, in the 
above expression for 07,, , . to make this correction 

^E(L<trciciu) 

(Hurford, D is always greater than unity, since an 
energy loss larger (smaller) than the mean energy loss 
early in the AE detector will result in a particle of 
lower (higher) than average energy and hence of higher 
(lower) than average specific Ionization. This particle 
wilt also tend to have a greater (smaller) than average 
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energy loss later in the detector. D is a function of 
L/R (the ratio of the AE detector thickness to the 
particle’s range) which increases from 1 at L/R = 0 to -^3 
for L/R'^0,‘15* The deceleration factor which we employ 
differs slightly from that used by Hurford since we 
include the velocity dependence of CPi£(LQn«iQtt7 obtaining 

this correct i on . 


So far the thickness, L, of the AE detector has not 
been specified. The effect which Landau fluctuations and 
energy measurement uncertainties can have on the mass 
resolution of the instrument can be minimized by an 
appropriate choice of this thickness. Figure 5 
illustrates this point. The long 45*^ line in figure 5a is 
a I ine of constant energy, 1550 MeV. The point along this 
line at which a given particle falls depends on the ratio, 
L/R, of the particle's path-length in the /lE detector to 
its range. The pairs of line segments intersecting the 
constant energy line are constant mass lines correspond i ng 
to the isotopes "^Be and as they would appear for 
various choices of the thickness L». The notable features 
of this graph are j) the separation between neighboring 
isotopes decreases as E' approaches E (or equivalently as 
L approaches 0) and 3^ for E ’ 0 the constant mass lines 
become parallel to the constant energy line while for 
E’^^E they become parallel to lines of constant AE. 


ORIGINAL PAGE IB 
OF POOR QUALITY 



3i 


FIGURE 5 

Example illustrating the effects of various 
choices of AE detector thickness, L, for the 
case of Be isotopes at a total kinetic energy 
of 1550 MeV ( ^Be range is HH-8 riini of Cs I ; '‘’Be 
range is Hi-1 of Cs I ) . 

a) Constant mass lines for the isotopes Be and 
'^Be for several values of the AE detector 

thickness. Separation between isotope tracks is 
greatest when Lisa large fraction of the 
part I c I e ' s range . 

b) Energy changes required to produce a change 
of i AMU in the calculated mass. The curves are 
I abe led as follows: 


AE 

-Change 

1 n 

AE 

with 

E» 

constant . 

E» 

-Change 

i n 

E ' 

w 1 th 

AE 

constant . 

Landau 

-Change 

1 n 

AE 

with 

AE+E' constant 


Also shown is the uncertainty in AE due to 
Landau fluctuations. 
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Figure 5^) shows, as a function of L, the amounts 
which AE and E* must separately be altered (as will be 
caused by errors in measuring these quantities) in order 
to alter the mass which we would calculate for the event 
by one mass unit. Also shown is the amount which />E must 
be varied with 4E+E* constant (as will be caused by Landau 
fluctuations) to produce a one mass unit change. In 
addition the figure shows the rms change in (and £') 
that will be caused by Landau fluctuations in the AE 
detector. It can be seen that the increasing separation 
of the mass tracks adequately compensates for the 
increasing Landau fluctuations as we increase the 
thickness of the 4E detector. 

In order to keep we utilize a stack of eight 
separate scintillation counters for energy measurements, 
The scintillator in which a particle stops is treated as 
the E’ detector and the energy losses in all the earlier 

V 

detectors are summed to form the /iE measurement. The 
thicknesses of the scintillators are chosen so that we 
obtain L/R^O-1 ' all hut the first few detectors. Table 
2 summarizes the thicknesses of the Csl(TI) scintillators. 

Figure b shows, for selected isotopes, the limiting 
mass resolution attainable with the HEIST instrument. It 
is assumed that for a particle stopping in Dn the energy 




TABLE 2 

HEIST Detector Thicknesses 


Detector 

' Norn 1 na I 

Th 1 ckness 

Sum Thru 

Th 1 s 

Maximum 




Detector 

L/R 


mm 

g/cm^ 

mm 

g/cm^ 


DO 

3 

1.4 

3 

1.4 

— 

Di 

3 

1.4 

b 

3-1 

0.50 

D2 

3 

1.4 

4 

4.1 

O.bl 

D3 

4 

1 8 

13 

5.4 

0.b4 

DH 

b 

3.1 

14 

8.b 

G.b8 

D5 

8 

3.b 

31 

13.3 

0.10 

Db 

13 

5.4 

34 

11. b 

0.b4 

1 

Q 

11 

1.1 ' 

5b 

35.3 

0.10 


LU 

i_n 



FIGURE b 


Contribution of Landau fluctuations to 
the uncertainty in measured masses for 
normally incident particles. All 
detectors prior to the detector in which 
the particle stops are summed to form 


the AE measurement. 



M, AE (LANDAU) 
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losses in DQ through Dn-i will be added to form the AE 
measurement. The sawtooth pattern is primarily due to the 
variation of the L/R ratio as discussed above. 

In addition to precise measurements of 4E and E’, it 
IS necessary to precisely measure the trajectories of the 
cosmic ray particles. Since a given fractional 
uncertainty in L, the particle's path-length in the /3E 
detector, results in the same fractional uncertainty in 
the calculated mass (see table 1) and since L is 
proportional to the secant of the particle's angle of 
incidence, B, it is essentia! to precisely measure d. 
Also, it is important to know the position at which each 
particle passed through the scintillator stack since 
corrections must be made for positional variations in the 
thickness and light collection efficiency of each 
scintillator. These corrections will be discussed in 
deta I I below. 


2. 2 Instrument — Overall Description 

Figure ”! shows a cross sectional view of the HEIST 
instrument. It consists of two basic parts. On top is a 
hodoscope consisting of eight multiwire proportional 
chambers (MWPC's) oriented so that alternate chambers 
measure x and y coordinates along a particle's trajectory. 



FIGURE 1 


Cross sectional view of the HEIST 
instrument. The cross section is taken 
along a diagonal of the proportional 
counter hodoscope. The light pipes and 
attached photomultipliers lie off of 
this plane and are shown projected onto 
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Beneath the hodoscope is a stack of nine thin, large 
diameter Csl(TI) crystal scintillators. The uppermost 
eight scintillators are used for making precise 
measurements of a particle's energy losses. Each of these 
scintillators is coupled to a 130 mm (5") photomultiplier 
tube by means of an adiabatic lucite light pipe. The 
ninth scintillator is directly viewed by a 38 
photomultiplier and is used to identify particles which 


d I d not 

stop 

in the 

upper 

eight 

sc 1 nt 1 1 1 ators . 

1 n 

addition. 

two 

plastic SCI, n til 

1 ators 

1 ocated 

be 1 ow 

1 

the 

cr ysta 1 

stack 

are used 

in anticoincidence to 

reduce 

the 


background of wide angle particles that can exit from the 
side of the crystal stack. Figures 8a and_b are schematic 
views from above and from one edge of the detector stack 
showing the relative spacing of the various detectors. 
The various sensors are identified as follows: Xj through 
XH are the MWPC's used for measuring the x-coord i nates of 
a particle's trajectory, Yi through YH are the MWPC's used 
for measuring the corresponding y-coord i nates , and DQ 
through D”| are the Csl(TI) scintillators used for making 
energy loss measurements. Table 3 lists the relative 
positions of the sensors shown in figure gb. 

Figure H illustrates, for particles with M=2Z, the 
intervals of charge and energy over which the HEIST 
instrument is capable of mass measurements. Two sets of 



FIGURE 8 


Scale drawings of particle detecting 
regions of the HEIST instrument. 

a) View from above. 

b) View taken through section A-A 
indicated in figure 8a- Note that this 
view 'is at H5° +o that shown in figure 
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TABL E 3 

Positioning of HEIST Sensors 

cm inches 


Bottom 

of 

hodoscope base plate 

0 

0 

Bottom 

of 

Y4 cathode 

0-435 

2-41b 

Bottom 

of 

X4 cathode 

i-853 

4-111 

Bottom 

of 

Y3 cathode 

M34 

b-451 

Bottom 

of 

X3 cathode 

3-bi8 

4-140 

Bottom 

o f 

Y2 cathode 

13-QB4 

33-081 

Bottom 

of 

X3 cathode 

13-401 

35-324 

Bottom 

o f 

Y| cathode 

14-184 

31-5b4 

Bottom 

of 

X| cathode 

15-bl8 

34-822 

Top of 

DO 

sc i nt I ! I ator 

-0-b48 

-l-b4b 

Top of 

Di 

sc I nt 1 1 1 ator 

-0-183 

-i-484 

Top of 

D2 

sc 1 n t i 1 1 ator 

-0-418 

-2-332 

Top of 

D3 

sc 1 nt i I 1 ator 

-l-04b 

-2-b51 

Top of 

DH 

sc 1 nt 1 1 1 ator 

-i-814 

-3-084 

Top of 

05 

sc f nt 1 1 1 ator 

-1-453 

-3-b41 

Top of 

Ob 

sc 1 nt 1 1 1 ator 

-1-110 

-4-44b 

Top of 

D1 

s c 1 nt 1 1 1 ator 

-2-245 

-5-102 

Top of 

D8 

s c 1 nt 1 1 1 ator 

-2-4ib 

-1-401 


* Reference level. 
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FIGURE S 


Intervals of charge and energy at the 
top of the earth's atmosphere which are 
accessible to the HEIST instrument for 
normally incident particles. Curves 
indicating the upper and lower limits of 
the energy intervals which can be 
detected are shown for atmospheric 
depths of 0, b and 1[) g/cm . 
Outlined regions correspond to flights 
at 1500 HV magnetic cutoff and b g/cm^ 
atmospheric depth (solid box) and at 
800 MV cutoff and 3 g/cm^ depth (dashed 
box) as discussed in the text. Energies 
correspond I ng to magnetic rigidities of 
800 and 1500 MV are indicated by arrows 
for particles with M/Z=P. 
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curves are included. The curves in each set are labeled 
according to the amount of overlying material (air plus 
material above the scintillator stack in the gondola). 
The lower set shows the energy which a normally incident 
particle must have at the top of the atmosphere in order 
to just penetrate to the OQ-Di boundary in the 
scintillator stack. The upper set ot curves is for 
particles which just penetrate to the D"|“D8 boundary. For 
normally incident particles of a given charge the energy 
interval above the atmosphere from which the HEIST 
instrument W’i I I accumulate stopping particles is simply 
the interval between the two curves corresponding to the 
appropriate amount of overlying material. We also 
indicate the energies which correspond to geomagnetic 
cutoff rigidities of 8QQ and 1500 MV In figure ^ we 
illustrate charge and energy regions viewed by the HEIST 
instrument in two cases. The solid box corresponds to a 
flight at a 1500 MV cutoff with b g/cm^ of residual 
atmospherev The relatively high magnetic cutoff in this 
case restricts observations to elements with Z 5: 5 and 
severely restricts the energy interval obtained for 
elements with Z^13. The dashed box corresponds to a 


* Such a flight was made on b June 1411 from Aberdeen, 

S.O. 



flight with a cutoff rigidity of gOQ MV at an atmospheric 
depth of 3 g/cm*^. In this case the cutoff has little 
effect and elements can be observed all the way down to 
the instrumental I imit of Z=3- It is clear from figure 3 
that it IS desirable to fly this instrument at northerly 
latitudes in order to obtain a geomagnetic cutoff 
<800 MV. 

In normal operation, a stopping particle must trigger 
at least one x and one y MWPC in both the upper and lower 
halves of the. hodoscope, must trigger two of the first 
three scintillators and not trigger either the bottom 
crystal or the plastic anticoincidence counter. Figure 10 
shows the integral geometrical factor (geometrical factor 
at angles less than a specified maximum angle) for 
accumulating stopping particles in two of the ranges of 
the instrument as a function of the particles* maximum 
accepted incidence angle. These geometrical factors are 


the result of a 

Monte 

Car 1 o 

ca 1 cu 1 at ion 

1 

(Sul 1 1 van , 

i^ni) 

and the error 

bars 

just 

reflect the 

finite number of 

trajectories used in 

the 

calculation. 

A I though 

the 

instrument is 

capab 1 

e of 

accepting events out to 

1 arge 

angles, the quality of 

the 

data will 

decrease at 

the 

larger angles 

due to 

increasing uncertainty in the 

angle 


of incidence and to the increased amount of atmosphere 
which these particles must traverse before reaching the 
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FIGURE 10 


Geometrical factor for collection of 
particles impinging on the detector 
stack at angles less than the indicated 
maximum angle. Geometrical factors are 
shown for particles stopping in 
detectors Dl and O']. The error bars are 
the result of the finite number of 
trajectories used in the Monte Carlo 
calculation. 
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I n str utnen t . 

Table ^ shows the number per day of cosmic ray 
particles of each of the elements from Li through Ni that 
are expected to stop in the instrument after traversing 
the overlying atmosphere without having undergone a 

nuclear interaction. Only events with angles of incidence 

o 

<30 from the axis of the instrument are included The 
values in table H are calculated using local 
interplanetary spectra based on the oxygen spectrum 
reported by Garcia-Munoz et al. (iSlla) and on the 
elemental abundance ratios compiled by Silberberg et al. 
from the work of a number of investigators. 

In table H rates are shown for geomagnetic cutoffs of 
800 and 1500 MV, typical of launch sites in southern 

Canada and the northern United States respectively, as 
well as for a variety of thicknesses of residual 
atmosphere. The effects of the residual atmosphere are 
twofold: first, ionization energy loss in the air causes 

a given energy interval at the instrument to map into a 
higher energy interval at the top of the atmosphere, 
thereby causing particles which had higher rigidity above 
the atmosphere to be observed, and second, nuclear 
interactions between cosmic ray particles and the air 
deplete the population of surviving cosmic rays while 



TABLE M 


HEIST Count Rates (events/day) vs. 
Atmospheric Depth and Geomagnetic Cutoff 



£ 1 ement 


800 MV 

cutoff 



1500 

MV cutoff 



Og/cm^ 

3g/cm^ 

bg/cm^ 

lOg/cm^ 

Qg/cm^ 

3g/cm*' 

bg/cm^ 

lOg/ 


L 1 


”412" 




”'Tr 



Be 

3bb 

331 

284 

230 

0.4 

38 

8b 

124 


B 

1200 

1010 

885 

b88 

134 

301 

424 

508 


C 

5330 

4b10 

3150 

2140 

5bl 

1240 

IbbO 

14b0 


N 

1510 

12b0 

440 

104 

33b 

480 

5bl 

b05 


0 

b2H0 

5020 

38b0 

2b10 

- 1400 

2350 

2510 

2b00 


F 

121 

48 

14 

50 

55 

bl 

b2 

50 


Ne 

1330 

1010 

143 

440 

548 

bOO 

bll 

440 


Na 

233 

no 

124 

80 

lib 

120 

118 

80 


Mg 

I1b0 

12b4 

844 

518 

84b 

813 

854 

518 


A! 

2b1 

184 

132 

84 

141 

141 

132 

84 


SI 

1380 

454 

b5b 

411 

135 

130 

b5b 

411 


P 

bb 

45 

30 

14 

34 

38 

30 

14 


S 

32b 

214 

141 

42 

188 

181 

141 

42 


C 1 

50 

33 

22' 

14 

31 

30 

22 

14 


Ar 

lb8 

no 

12 

45 

102 

41 

12 

45 


K 

% 

b2 

40 

25 

b2 

58 

40 

25 


Ca 

2b5 

lb8 

no 

b8 

lb4 

158 

no 

. b8 


Sc 

44 

28 

18 

11 

31 

28 

18 

11 


T 1 

210 

132 

84 

51 

155 

132 

84 

51 


V 

40 

5b 

35 

21 

b8 

5b 

35 

21 

o o 

Cr 

185 

113 

12 

43 

134 

113 

12 

41 


Mn 

111 

11 

45 

21 

40 

11 

45 

21 

si 

Fe 

1440 

8b2 

545 

324 

1103 

8b2 

545 

324 

Co 

b.l 

3-b 

2-3 

1-3 

4.8 

3-b 

2.3 

1.3 

€s 

N 1 

b3 

31 

23 

14 

48 

31 

23 

14 

tot a 1 

23300 

18400 

14000 

4800 

. 1b00 

8800 

4000 

8300 

si 

events/ sec 

.21 

.21 

.lb 

'.11 

.088 

.10 

• 10 

• 10 
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producing atmospheric secondaries. The fraction " of the 
particles which have undergone nuclear interactions at 
various depths in the atmosphere can be seen by comparing 
the rates listed for a 800 MV cutoff in Table H for these 
depths with the rate with no residual atmosphere. 
Calculations by Silberberg and Tsao (1S11) indicate that 
the isotopic composition of most elements does not vary 
greatly with atmospheric depth. However, It is clearly 
desirable to minimize the fraction of atmospheric 
secondaries collected in order to reduce the uncertainty 
involved in extrapolating abundance ratios to 
interplanetary space. 

The total thickness of the crystal stack, 
25-3 g/cm^ , IS not negl igible compared to the 5b g/cm** 
nuclear interaction length of Fe in Cs I . In the longer 
ranges approximately 1/3 of the Incident iron particles 
Will undergo a nuclear interaction and become unusable for 
mass deterifl I na t i on . This limit excludes the possibility 
of extension of the ^E-E’ technique to measurement of 
cosmic ray isotopes at energies greater than a few 
GeV/n uc . 
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E-3 Tr,ajectory Measurement 

The hodoscope consists of eight identical multiwire 
proportional chambers. Figure H shows a detailed view of 
a single chamber and table 5 summarizes the mechanical 
characteristics of the chambers. The sensitive area of 
the chamber is H8cm x 48cm (8300 cm^ area). The frames on 
which the proportional counter planes are constructed are 
machined out of G-IQ epoxy-glass material. Printed 
circuit boards with copper pads to which the chamber wires 
can be soldered are inlaid in the frames. The anode plane 
IS composed of .80 y^m diameter stainless steel wires spaced 
4-8 mm apart. There are two cathode planes, each spaced 
4 mm away from the anode plane. The upper cathode 
consists of a sheet of b*4 thick aluminized mylar and 
serves only to establish the appropriate voltage 
difference across the 4 mm cat hode-to-anode gap. The 
lower cathode is composed of 85 /•m diameter stainless 
steel wires spaced ?*1 mm apart. This plane is used both 
to obtain the required potential difference and to sense 
(in one dimension) the location at which the avalanche 
takes place. The cathode wires are oriented perpendicular 
to the anode wires in order to allow position measurements 
with finer resolution than the spacing between wires, as 
will be discussed below. The bottom plane in each MWPC is 
another sheet of b ■ 4 aluminized mylar. This plane. 
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onal counter. Thicknesses are 
4x over the indicated scale for 
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MWPC Mechanical Characteristics 


frame mater i a ( 


G-10 ( epoxy~g I ass ) 


active area 


4S cm X 48 cm 


cathode-to-anode spaci 

1 ng 

4 mm 


cathode-to-suppressor 

spacing 

8.4 mfn 


anode wire spacing 


4*23 fnm 


cathode wire spacing 


5-12 nfim 


anode wires 


20 diameter stainless 

stee 1 

cathode wires 


25 diameter stainless 

stee 1 

upper cathode and suppressor 

b-4y«nfi mylar, aluminized 

on both sides 
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which we call a suppressor, is operated at approximately 

3/4 of the anode voltage and serves to collect charge 

✓ 

deposited in the region below the cathode wires. We 
employ a cathode coupled delay line readout of the 
Perez-Mendez design (Grove et al., i4"10» 1412* 1413)* 

The operation of multiwire proportional chambers has 
been reviewed by Charpak (1410)- Electrons which are 
detached from atoms of the chamber gas by the passage of a 
charged particle through the sensitive volume of the 
chamber rapidly drift toward the nearest anode wire. In 
the region of large electric field near the anode wire 
they avalanche. The resulting secondary electrons are 
collected on the anode wire, and as the remaining sheath 
of positive Ions drifts slowly (M50 nsec) through the 
region of large field near the anode wire an image charge 
IS induced on the cathode wires. The centroid of this 
image charge .distribution in the direction perpendicular 
to the cathode wires represents the location at which the 
cosmic ray particle traversed the chamber since the 
primary electrons have not drifted in this direction. 

The image charge is distributed on the cathode wires 
approximately according to 

^ _J 
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where 1^ is the ca thode-to-anode spacing, is the 

location of the n"^ cathode wire and x is the x-coordinate 
of the cosmic ray trajectory (Lacy and Lindsey, iS13)* 
This expression is easily derived by integrating the image 
charge distribution on a conducting plane (the cathode 
plane) due to a point charge (the avalanche) located, a 
distance 1^ above the plane along parallel strips of equal 
width (the wires). We have measured this distribution in 
a prototype of our proportional chambers by pulse height 
analyzing signals on a single cathode wire produced by a 
collimated source of 5*3 keV x-rays ( Fe). The distance 
of the source from the wire being analyzed was varied in 
precisely known steps. The resulting variation of pulse 
height with source position cap be seen in figure IJ. In 
addition, the predictions of equation b are shown. We 
have used the independently determined value of 
but have adjusted the height and center 
position of the curve to fit the data. 

For our flight chambers 1^ is 4 mm and the spacing 
between adjacent cathode wires is 2-i so i’he 

total image charge is distributed over cathode wires. 

The centroid of this distribution can be interpolated to a 
fraction of the wire spacing. The cathode signal is 
capacitivety coupled onto an electromagnetic delay line, 
resulting in a pair of pulses travelling in opposite 
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Distribution of image charge on MWPC 
cathode wires due to avalanche at the 
anode The data were accumulated using a 
collimated ^®Fe source of 5.S keV 
x-rays Error bars represent 
uncertainties in determining the mean of 
the accumulated distributions of pulse 
heights. The sol id curve shows the 
expected shape of this distribution, as 
discussed in the text. 



PULSE HEIGHT ON SINGLE CATHODE WIRE 
(ARBITRARY UNITS) 
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directions down the line. Measurement of the difference 
of the arrival times of these two pulses allows one to 
determine one coordinate of the cosmic ray's trajectory. 

The delay lines which we use are modeled after those 
developed by the Perez-Mendez group (Grove et al., iS13» 
Perez-Mendez and Parker, 19145 • Table b lists the 
characteristics of our delay lines. It should be noted 
that the delay and attenuation measurements were obtained 
while the delay line was clamped to the printed circuit 
boards of the chamber cathode, thus they include any 
effects due, for example, to added capacitance to ground. 
Figure 13 shows pulses obtained at one end .of the delay 
line from identical signals injected at three positions 
along the line-*-near each of the ends and near the center. 
It is clear from the picture that not only does the pulse 
suffer attenuation by a factor It traverses the 
entire length of the line, but i+s shape is noticeably 
altered as well. The timing circuitry requires that the 
signals received at the two ends of the delay line be of 
approximately the same shape and amplitude. The variation 
shown in figure 13 caused excessive variation of the time 
difference obtained by these circuits as > the signal 
amplitude was varied at a fixed position along the delay 


line. 



Tabl e b 

Delay Line Characi" eristics 


Mech a n i ca 1 

f ength 
w I dth 
th I ckness 
core mater i a I 
ground plane 
compensation material 

winding 


53 cm 
3.2 cm 
0.33 cm 

G-1 (fiber glass with silicone binder) 
copper strips on kapton backing 

aluminum strips on mylar backing, o- 

strips at 45*’ to axis of I ine 

8? turns/cm (808/inch) of 38 ^WG magnet wire 


Electrical 


de 1 ay 

characteristic impedance 
attenuation (for MWPC pul 
DC res 1 stance 


D O 

tog 

s >■ 

^£3 


ses ) 


15 nsec/cm 

i50on 

factor of 8 P®f' 3H cm 
lIXL/cm 
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F I GURE JJ 


Pulse shapes received at one end of 
53 cm long delay line for identical 
pulses injected near the receiving end 
(a), near the center (b) and near the 
opposite end (c). Tick marks on the 
horizontal axes are separated by 
100 nsec and those on the vertical axes 
are separated by 4 mV. 
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To reduce the magnitude of this effect a system of 
taps was devised. The delay line is divided into thirds 
as shown in figure 14. Signals from one end of the delay 
line and from the delay line tap located nearest to the 
opposite end of the line are connected together at the 
input to a preamplifier located midway between these 
points. Similarly signals from the other end and the 
other tap are connected at the input of a second preamp. 

In obtaining the time differences for an event the timing 
circuitry ignores all but the first pulse from each 
preamp, so the position measurement is unaffected by the 
eventual arrival of the delay line pulse at the second 
input to the preamp. Also, reflections from the ends of 
the delay line (the preamps have input impedances much 
less than the characteristic impedance of the line) will 
not affect the measurement. 

* 

With this arrangement of taps the maximum time 
difference corresponds to 1/3 of the entire delay line and 
the degree of attenuation and distortion of the pulses is 
reduced. However, an ambiguity has been i ntroduced--any 
given time difference will correspond to three different 
positions in the chamber. In order to resolve this 
ambiguity a system of charge sensitive amp I ifiers and 
comparators was introduced. A fixed time after the 
passage of a charged particle through the instrument the 

ORIGINAL PAGE .. 
OE'POOR QUALm 
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F I SURE IH 


Block diagram of position sensing 
circuitry. 




a~ 
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outputs of these sector preamps are examined and the 
number of the sector having the largest signal is encoded 
and inserted into the data stream for that event. 

Figure iH is a block diagram of the position sensing 
circuitry and figure 15 shows how a digital time 
difference is obtained from the two delay line pulses. 
The output of a delay I i ne preamp goes into a "high level 
amplifier" which differentiates and amplifies the pulse. 
The output of this amplifier goes two places; first, to a 
discriminator which will produce a logic pulse at the time 
of the zero crossing of the differentiated pulse if the 
pulse amplitude is sufficiently large, and second, to 
another stage of amplification followed by an identical 
discriminator. The two discriminator pulses are "ORed" to 
produce the pulse which defines one end of the time 
interval. This dual amplifier and discriminator 
arrangement is required to handle the large dynamic range 
of signals expected from the hodoscope, ojilOsi* 

When the pulse is large enough to trigger the high 
level discriminator, the timing will be done using this 
pulse since propagation delays through the low level 
amplifier ensure that the high level pulse will arrive 
first. The discriminator pulse initiates the charging of 
a capacitor by the time stretcher circuit. When the next 



1i 


F I GUR E J^5 


Example illustrating the processing of 
delay line pulses to obtain a digital 
t ( me d i f f erence . 



ORIGINAL PAGE IS 
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;jj¥injuuiiiRjirijmiinjiJiiiii^^ running 

-^(^I5nsp.n t 33,33 MHz CLOCK 




ISnsec 

DELAY LINE PULSE A 

DIFFERENTIATED PULSE 

DISCRIMINATOR LATCH 
(produced at zero crossing of differentiated 
pulse) 

A TIME 
STRETCHER 
RAMP 





15 PULSES 
GATED TO 
A- SCALER 


yiv 


'810 nsec 


A = I5 
B = 7 
C = 27 

TIME DIFFERENCE = 30C-fA-B 

= 8l8nsec 


nsec 

DELAY LINE PULSE B 
DIFFERENTIATED PULSE 
ZERO CROSSING LATCH 




nnjuTJTJTJiAjmn^ 


27 PULSES 
GATED TO C SCALER 


B TIME STRETCHER 
RAMP 

7 PULSES GATED 
TO B SCALER 
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pulse from the 33.33 MHz clock arrives the time stretcher 
starts to discharge this capacitor at l/pS of the charging 
rate. During the discharge time clock pulses are gated to 
the "A-scaler". The final count in this scaler represents 


the time 

i nterva I 

(measured in 

nanoseconds) between 

the 

ar r I va I 

of the 

pulse from 

the 

delay line and 

the 

occurrence of the 

next c I ock pulse. 

The other delay 

line 


pulse IS treated in an identical manner resulting in a 
second vernier measurement in the "S-scaler". Finally, 
the ”C-scaIer" counts the number of clock pulses between 
the arrivals of the two delay line pulses. This C-scaler 
count gives the gross time interval (in units of 30 nsec), 
and the combination 3QC+A-8 is the total time difference 
in nanosecQnds. 

We have found that in ordet to obtain position 
resolution which is not limited by the spacing of the 
cathode wires it is essential that the charge deposited 
outside of the active volume of the chamber, beyond the 
plane of cathode wires, not be al lowed to drift to the 
anode. The suppressor plane was introduced to achieve 
this result. If this charge were permitted to drift to 
the anode, it would distort our measurements of the 
position of the cosmic ray's track since the electric 
field I ines which the electrons f o I low in this region have 
a component in the direction which we are trying to 
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measure. This effect is illustrated In figure lb* A 

ss 

prototype chamber was illuminated by an uncoil imated Fe 
source and distributions of the resulting positions were 
accumuiated. The data in figure iba were accumulated with 
the suppressor grounded. The resulting distribution 
consists of a number of peaks produced by the effect 
described above added to a smooth distribution obtained 
from x-rays which interacted in the normally active volume 
of the chamber. Data in figure Ibb were taken with the 
suppressor biased at 1500 volts. The peaks have been 
eliminated since the suppressor can now collect the charge 
deposited in the volume between the cathode plane and the 
suppressor. in both distributions the smooth large scale 
structure is due to non-uniform illumination by the x-ray 
source. When the chamber is illuminated with charged 
particles the effect of the suppressor is less dramqtic, 
since in this case each event deposits charge in both the 
active region and in the region between cathode and 
suppressor planes. When the suppressor is not biased the 
measured position wi-ll be an average of the positions at 
which the charge from the normally active region and the 
charge from the region between the lower cathode and the 
suppressor arrive at the anode of the chamber, weighted by 
the magnitudes of the charges. 
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F I GU RE lb 

£5 

MWPC response to uncol I imated Fe x-ray 
source. The spacing between adjacent 
cathode wired is indicated. 

a) No bias applied to suppressor. 

b) Suppressor biased at 1500 vol+s. 



400 



POSITION IN MWPC (mm), ARBITRARY ZERO 
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We use a mixture of 10? Ar plus 30? as the 
proportional chamber gas. We have measured the gain 
characteristics of our chambers when stimulated with 
various sources of radiation. Figure H shows the charge 
collected on the chamber anode as a function of the anode 
voltage for each of these radiations. Where error bars 
are included they Indicate the range of charges collected 
when a number of different proportional counters were 
tested.' These gain differences are believed to be due to 
small variations in the chamber geometry (which will alter 
the chamber capacitance and thereby alter the static 
charge on the anode and the chamber gain) and possibly in 
part to a small degree of contamination of the chamber 
gas. 


Our estimates of the energy deposited in the 
i*3H mg/cm^ active thickness of the chamber are indicated 
in the figure. At any given value of the anode voltage, 
the measured charges do not increase in proportion to the 
energy deposited. The dashed lines in the figure indicate 
the magnitude of the static charge on sections of the 
anode wire Q.5 mm and l.Q mm in length. Doolittle et al. 
^1313^ estimate that the effective length along the anode 
wire of the distribution of charge produced in the 
avalanche will lie within these limits for a chamber 
geometry similar to ours. When the charge in the 
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FIGUR E J^l 

Charge collected on MWPC anode when 
exposed to various sources of radiation 
as a function of anode voltage. The 
error bars, where present, indicate the 
range of measurements obtained from a 
number of chambers that were tested. 
Also indicated are the charge levels 
(indicated by dashed lines) at which 
space charge effects are expected to be 
significant, as discussed in the text. 



CHARGE COLLECTED ON ANODE (pC) 
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avalanche is comparable to this static charge, a reduction 
in the gain of the chamber is expected since the charge in 
the electron cloud reduces the electric field seen by i ate 
arriving electrons. 

Ideally, the position measurements made using the 
techniques described above should be independent of the 
amplitude of the pulse being measured. However we observe 
a residual peak-to-peak time variation of <^i5”25 nsec 
when we st i mu I ate. the chamb’ers with pulses of various 
amplitudes at a fixed position. The main dependence on 
pulse amplitude of the technique which we use to measure 
position is caused by failure of the analog circuits at 
the two ends of the delay line to precisely track one 
another due to attenuation of the delay line pulses. 

Figure 18 shows the transfer function (time vs. 
position) obtained in one of the hodoscope planes. These 
data were obtained by accumulating distributions of times 
while exposing the chamber to a beam of collimated alpha 
particles obtained from the decay of Bi (a short lived 

A AiS 

descendant of Th). This source was mounted on a 
movable stage driven by a precision lead screw. After the 
desired number of events (typically 4000~5000) had been 
accumulated at a given source position, the source was 
moved a selected distance and the next distribution was 
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FIGURE JJ^ 

MWPC transfer function (time vs. 
position) measured using collimated 
c<-source. The errors in measuring both 
the position and time difference are too 
small to be displayed on this scale. 



MEAN TIME DIFFERENC 



SOURCE POSITION (cm), ARBITRARY ZERO 
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accumulated. The figure shows that there is a sizeable 
deviation from linearity near each of the taps 

Distortion of the transfer function is to be expected at 
positions located closer to a tap than the width of a 
typical delay line pulse, (''i50~800 nsec or A.-10-13 mm). 
At these places pulses travelling in both directions 
arrive at the most distant delay line preamp close enough 
in time to overlap and cause the shaping amplifier to 
slightly alter the time of the zero-crossing of the 
differentiated delay line pulse. 

Figure 18 shows a more detailed view of the response 
of one pair of hodoscope planes over a 10 mm interval near 
the center of the chamber. These data were accumulated at 
source positions 0-2 mm apart Figure 18a is a plot of 
time vs. position, figure 18b shows the standard 
deviation of the time distributions obtained, and figure 
18c shows the deviations of the means of the distributions 
from a least squares straight I ine fit. The fact that the 
deviations from a straight I i ne are smal I, at least over 
distances as large as 1 cm, shows that a relatively small 
number of cal ibration points are required to obtain the 
transfer function of each hodoscope plane to the accuracy 
that we requ i re . 
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FIGURE iS 


MWPC 

response 

over a 

1 1 m 1 ted sect i on 

of 

the 

chamber 

( i cm ) , 

measured w i th 

6 

CO 1 I 

1 mated o<- 

source . 




a) Mean time (yu) vs. position transfer 
f unct ion. 

b) Standard deviation ((T) of the 
measured position distributions. 

c) Deviation (A) of transfer function 
from best-fit straight line. 



o-(nsec) 
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SOURCE POSITION (mm), ARBITRARY ZERO 
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Figure 20 shows the MWPC resolution (rms) as function 
of signal amplitude. The points indicated by dots and 

triangles were obtained by using a pulser signal with a 
150 nsec rise-time Injected onto one of the cathode wires. 
The shape of this curve is the result of random noise in 
the input stages of the delay line amplifiers. Also 
included on the plot are points obtained by using a 
collimated ^''’Po alpha source and varying the chamber high 
voltage in order to vary the pulse amplitude prociuced. 
The correspondence between pulser amplitude and anode 
charge was obtained by observing pulses from both the 
pulser a'nd from charged particles at the output of the 
delay line preamps. The rms widths of the 0( 
distributions track those of the pulser well at low 
amplitudes, but then level off at nsec 

{FWHM^9-4 nsec). This is predominantly the result of the 
finite size of the collimator being used. 

In figure 2i we show typical time distributions 

acquired using various sources of radiation. A collimated 

source of alpha particles (8*185 kinetic energy) from 

the decay of Bt was used to produce the distribution in 

figure 2i9* Figures 2ib and c show histograms of time 

differences of x-coord i nates measured in the X| and X2 

proportional counters when the hodoscope was exposed to 

'lo 

parallel beams of fully stripped Ar ions at energies of 
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FIGURE 20 


Dependence of the rms position 
resolution in a single MWPC on pulse 
amplitude. Measurements are shown for 
the chamber stimulated with .an 
electronic pulser and with a collimated 
source of ©(-particles. The insert shows 
the collimator geometry used. 
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FIGURE 


Time distributions illustrating MWPC 
resolution. 

a) Chamber stimulated by a collimated 
source of <X- particles. 

b) Chamber stimulated by a parallel beam 
of normally incident, fully stripped 

hr ions at 810 MeV/nqc. The abscissa 
Is the difference between position 
measurements in the Xl and X8 
proportional counters. The areas of the 
hodoscope from which events are selected 
is discussed in the text. 

4o 

c) Same as (b) for J.55 MeV/nuc hr. 












810 MeV/nuc and J.55 MeV/nuc respectively. Events selected 

for inclusion in these plots lie within a band extending 

2000 nsec {'^13 cm) in the y direction and 33 nsec (~2 

’in. the X direction. By using the difference Xl-X2» we 

e! iminate the mean variation of the measured time over 

this 33 nsec x interval. Since both XI and X2 are subject 

to random errors and these errors are independent and are 

assumed to come from distributions of the same standard 

deviat.ion, we must divide the standard deviation obtained 

from the distribution of Xj-X^ by 4?* ’ order to obtain 

the resolution of a single hodoscope plane. The 

to 

resolution obtained using energetic Ar ions is 

significantly larger than that obtained using alpha 
particles (where the resolution obtained is limited by the 
size of the collimator employed). This may be related to 
the production of knock-on electrons in the former case. 
However, this possibility has not been thoroughly 
investigated. 

The traj'ectory of a particle passing through the 
hodoscope can be determined from two measurements of its x 
position and two measurements of its y position along its 
track. For each event we make four measurements of x and 
four of y. These redundant measurements permit a 
consistency check of the hodoscope data. In addition, 
doubling the number of position measurements improves the 



overall position resolution by a factor of approximately 
V? . In appendix B the uncertainties in sec0 and in the 
position at which the particle's track passed through the 
various scintillators are derived from the hodoscope 
resolution and the geometrical configuration of the 
instrument. For the hodoscope geometry used in this 
experiment the relative uncertainty In our determ i n at I on 
of sec 0 IS found to be 

..,11?.^ zz O.OO/^ ' 

where ^ is the particle's angle of incidence and Is 
the position resolution in a single proportional counter 
in millimeters. This contribution to our thickness 
uncertainty is shown in figure 2E* 

Also shown In figure 22 Is the contribution to our 
thickness uncertainty due to the crystal thickness 
variations (to be discussed later). This contribution 
depends on the precision with which we can determine the 
position at which the particle passed through each of the 
sc 1 n 1 1 I I ators . In appendix B we show that the uncertainty 
In the X and y positions at which a particle’s track 
Intersects a given level in the scintillator stack varies 
between 0*81 and l.QS times . 
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FIGURE 22 


Relative uncertainty (rms) in AE 
detector thickness, L, along with 
contributions due to angular uncertainty 
and to thickness mapping uncertainty. 
These uncertainties are plotted as a 
function of the particle's angle of 
incidence, 





L (thickness) 





2*4 Energy Measurements 


A detailed view of the scintillator stack is shown in 
figure P3. Each scintillator is a disk of Csl(TI) 
(Harshaw Chemical Company, Solon, Ohio) eleven inches 
^Bl-4 cm) in diameter. The thicknesses vary from 3 
^•35 g/cm*‘) to 13 inm (l.bl g/cm^) . The top and bottom 
faces of each disk have Harshaw’s "standard crystal 
polish" while the circumferences are unpolished. Each 
crystal is connected to a lucite light pipe. The light 
pipe is coupled to the crystal along 1/2 of the crystal's 
circumference through a RTV~b02 transparent silicone 
rubber gasket. The opposite edge of the crystal is in 
contact with a thin strip of Millipore fitter paper (type 
HAWP) which acts as an efficient diffusely reflecting 
surface at the wave lengths of the Csl(TI) emission. This 
filter paper is backed by another RTV-bOE gasket which 
acts as a deformable cushion attached to a PVC plastic 
frame. A pair of screws pull this assembly together and 
force the light pipe into optical contact with the 
scintillator. The light pipe is divided into six lucite 
strips which are bent so as to adiabaticly transport light 
to the face of a 130 mm (5") photomultipl i er tube (EMI 

S530R). 
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FIGURE 23 


Mechanical assembly of the scintillator 
stack. 

a) Coupling between CsKTI) scintillator 
and lucite light pipe. 

b) Material used to separate adjacent 
crystals (not to scale). 




Successive scintillators are separated by two sheets 
of Millipore and a single sheet of aluminized mylar. The 
total amount of material between crystals is approximately 
1 mg/cm^. The aluminized mylar is used to ensure that 
light produced in one crystal will not be transmitted into 
an adjacent crystal. The Millipore sheets in contact ‘with 
each of the crystal faces improve the transport of light 
to the light pipes. Tests indicate that diffuse 
reflection by the Millipore dominates internal reflection 
at the crysta|~air interface, presumably due to inadequate 
polishing of the crystal surfaces. On some of the 
crystals an additional "compensation piece" of aluminized 
mylar was added over a fraction of the crystal surface 
between the crystal and the Millipore. This piece serves 
to reduce the light collection efficiency in areas where 
it otherwise would be excessively large. The result is 
smaller light collection gradients. 

It IS necessary to obtain thickness maps of each of 
the scintillators. Low resolution maps were made for each 
of the crystals by looking at the variation of the 
attenuation of a beam of collimated, low energy ^-rays. 
The crystal to be measured was positioned betv/een the 
^-ray source ( with = keV or ^^^Ba with 
Ejf = 35b keV) and a small Nal(TI) scintillator (?.5H cm 
diameter x 2-5H cm thick) mounted on a pm tube. The 



counting rate in the photopeak seen in the NalCTI) 
scintillator was measured with the Csl(TI) crystal in a 
number of positions. The difference in thickness between 
two places on the crystal can be obtained from 

^ ^ Jin (jf) 

where 4z IS the thickness difference, is the ^-ray 

absorption coefficient of Csl(TI) at the energy of the 
2^-ray being used, and Nj, and N are respectively the number 
of counts obtained in a fixed time interval at a reference 
position (the- center of the disk) and at the position of 
interest. The number of counts required in order to 
obtain an uncertainty of in the thickness difference 

measurement i s 

N= 

where we have used the fact that for the thickness 

differences which we are considering. In order to obtain 
thickness differences with uncertainties less than Q.i^ of 
the crystal thickness it was necessary to obtain between 
1x10 and 8x10 counts at each position, depending on 
which crystal was being mapped and which of the ^-ray 
sources was being used. Each thickness measurement was an 
average over an area "8 to 3 cm i n diameter. rt'ff TS 
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One of the resulting maps Is shown In figure 24* The 
maximum thickness gradients obtained over the area which 
we were able to map ranged between 1.2 and 2 • 4 on 
all crystals except 05. The 05 crystal exhibited 
thickness variations as large as b ^m/mm. If we estimate 
the thickness gradients in the various combinations of 
crystals to be used for the measurement (that is, in 
the combination OQ through Dn-l for the case where a 
particle stops in On) by simply summing the maximutri 
thickness gradients in each of these crystals, we find 
that this overall gradient is less than 0-i^ of the AE 
detector thickness per mm in all cases. This gradient 
combined with our position resolution yields a relative 
thickness uncertainty of 


< o.oof 

~L 

Here we have assumed that uncertainties in the thickness 
at any given point in the crystal can be calibrated with 
an uncertainty <0-03? so that the uncertainty in 
obtaining thicknesses will only depend on the precision 
with which we can determine the point at which a particle 
passed through the crystal and on the thickness gradients. 
This condition should be obtainable from a heavy ion 
calibration with reasonable counting statistics. In 


figure 23 i+ can be seen that the thickness uncertainty is 
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Thickness contour map 
scintillator. Thicknesses 
differences from thickness at 


of Dl 
shown are 
the center 


of the crystal. The nominal thickness 
of this crystal is 3 "i"!* Thickness 
gradients are also indicated. 
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comparable to the uncertainty In L due to the finite 
angular resolution of the hodoscope. 

Figure ?5 shows the nominal response of the 

scintillators to various normally incident particles 
stopping in D"| . The maximum energy lo.ss of interest in D*1 
IS 20 GeV for normally incident particles and 22 for 

particles incident at 30*- A+ low energies we shou.l d be 
able to detect particles which have penetrated only a 
small distance into Dl . Failure to detect such particles 
would cause them to be treated like particles which 
stopped in Db and would introduce an error into our 
measurement of E*. If a particle were to penetrate 10 
into D3 and not be detected, a 0-b^ error in our estimate 
of the energy loss in Db (E') would result. This amount 
of material is a small fraction of the total thickness of 
the E* detector and therefore few particles will stop in 
this layer, so it is not unreasonable to exclude from 
analysis those particles which stop sufficiently near the 
boundary of two detectors that it Is unclear in which 
detector they stopped. If we do exclude such events, the 
sma I lest signal of interest will be that produced by our 
^-ray calibration sources, to be discussed later. The 
calibration source mounted on the Dl pm tube produces 
pulses approximately equivalent to an energy loss of 
l.b in the Dl crystal. The overall dynamic range 



FIGURE P5 


Nominal AE vs. E’ tracks for 
incident particles stopping 


norma I I y 
in Dl are 


shown by solid lines. Lines of constant 
particle range are also shown (dotted). 
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required in Dl (assuming no saturation of the 
scintillation at large energy losses) Is 

Figure 2b shows the circuitry used to analyze the 
photomultiplier signals. Dynodes 8 through J.J, of the 
photomultipliers are connected to the anode since for the 
large signals produced by heavy ions a large tube gain is 
not required. Signals from dynodes 3 connected 

f 

to separate ampi ifier chains. This arrangement Is 
necessitated by the large range of energies deposited in 
the scintillators by the cosmic rays in which we are 
interested. The photomultiplier gain up to the stage 
being, analyzed must be large enough so that electron.ic 
noise wi I I not produce a large uncertainty in our energy 
measurements, even for the smallest signals of interest. 
However, we do not want to analyze signals from a dynode 
at which extreme space charge saturation is occurring. 
Due to the long decay time of the Cs!(TI) scintillation 
{ r-»i yusec) relatively large amounts of charge can be 
delivered to the dynode without saturation occurring. 
However, the onset of saturation occurs near the upper end 
of the range of pulse heights which we analyze op each 
dynode . 
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FIGURE 2b 


Block diagram of pulse height analyzer 
circuitry. 
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Figure 21 illustrates the saturation characteristics 
of dynodes 3 and "| in one of our photomultipliers. In 
this figure the ratio of the charge, Q, collected at the 
dynode of interest to the intensity of the light pulse, !, 
incident on the cathode of the tube Is plotted as a 
function of 0 . The tube was biased with H50 volts between 
the cathode and dynode 1 (for focusing purposes) and with 
80 volts between each adjacent pair of dynodes up through 
the dy'j-dyS pair. Dynodes 8 through H were connected to 
the anode. This biasing arrangement yielded a ratio of 
iOS between the charges collected at dy"] and dy3 for a 
light input small enough so that dy1 exhibited no 
saturation. The cathode of the tube was illuminated by a 
light emitting diode (LED) pulsed in such a way as to 
produce an approximately exponentially decaying light 
pulse with a lyAsec time constant. The LED was placed at 
the focal point of a lens situated between the LED and the 
face of the photomultiplier. The lens served to produce 
an approximately parallel beam of light incident on the pm 
tube face. Between the lens and the tube face were 
inserted various combinations of neutral density filters 
for varying the light intensity. Each combination of 
filters used was calibrated using the pm tube response in 
the region where no space charge effects can be observed. 
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F I GURE n 


Saturation characteristics of EMi S53QR 
pm tube signals. The bias network used 
IS discussed in the text. Uncertainties 
in the plotted points are dominated by 
systematic uncertainties in determining 
the light intensity incident on the pm 
tube cathode. 

a) Saturation of signals measured on 
dynode I* Open circles indicate data 
taken after set-up was , modified to 
permit measurements at I arger pulse 
heights. 

b) Saturation of signals measured on 
dynode 3- 
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As can be seen in figure 21 the onset of space charge 
effects occurs at dynode 1 for a charge /u10 pC and at 
dynode 3 tor a charge n-'lO pC. This difference Is 
probably the result of the large space charge density at 
later dynodes when dynode 3 col lecting charges "^10 pC 
(when dynode 3 is collecting Ifl pC, dynode 1 will be 
collecting ^^550 pC). It should be noted that we are able 
to operate our pm tubes with a certain degree of 
saturation without significantly degrading the performance 
of the instrument. Since we are using logarithmic pulse 
height analyzers to compress the signals, additional 
compression due to pm tube saturation can be tolerated if 
Its effect is small compared to the logarithmic 
compression by the analyzers. Note, however, that pm tube 
saturation will complicate the i nter pretat i on of our pulse 
height measurements. in particular, a gain change (due, 
for example, to a change in the pm tube’s temperature) 
will produce a different factor change in the observed 
signal size for pulses which cause pm tube saturation than 
for those which do not. We estimate that this effect will 
alter the pm tube temperature coefficient by for the 
largest signals of interest. 

The analog circuits used to analyze the pm tube 
signals have a dynamic range of approximately inO^l* 
Also, the voltage distributions on the dynodes yields a 
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ratio ^lOQsl between the signals obtained from dynodes "I 
and 3- Ori each dynode signals In the range 1 pC to IDO pC 
are analyzed, so that an an overall dynamic range of 
/^lO ;1 IS obtained. Logic circuitry determines which of 
the two dynode signals will be pulse height analyzed based 
on whether the dynode 3 signal is large enough to trigger 
the d I scr I m I nator . The pulse height analyzers have a 
logarithmic transfer function with an average channel 

widthof0»133^- 

V 

The process of detecting and encoding the energy 
losses in the scintillators involves the following steps; 
1) production of light as ion-electron pairs produced by 
tho passage of a charged particle are recombined, 
trans.port of light from the point of production the 

light pipe, 3^ transport of light down the light pipe to 
the face of the pm tube, production of photoelectrons 
as the light impinges on the photocathode, 5^ gain by 
means of secondary emission at each of the dynodes, 
amplification and shaping in the analog stages of the 
pulse height analyzer, 3^ digitization in the 
ana I og-to-d i g I ta I converter. Each of these stages 
introduces additional uncertainty into the measurement. 
The major sources of uncertainty are described in the 
following sections. Figure shows the energy dependence 
of each of^ these sources of uncertainty in the energy 
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Energy 

measurement 

uncerta i nty 

1 n a 

single 

detector vs 

. energy deposited. 

Various 

contributions to the 

overa 1 I 

energy 

measurement 

uncertainty 

are also 

shown ^ 

as described 

in the text . 
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measurement. In addition the overall energy measurement 
uncertainty is shown. 


The transport of photons to the photocathode and the 
product'ion of electrons by photoelectric absorption and 
then by secondary emission constitute a chain- of 

statistical processes. The compounding of the statistical 
errors in such a chain Is treated. in detail by 

Br e 1 tenber ger (iS55)* The overall statistical uncertainty 
IS dominated by the stage where the signal is being 
transported by the fewest carr i er s--that is by the 
production of photoelectrons at the cathode. The 
fractional uncertainty in the signal is approximately 
where N is the number of carriers ( photoe I ectrons ) 
at this stage. We obtain a photoelectron yield of 
approximately 150 per MeV deposited in our scintM lators. 
This estimate of the yield is uncertain by a factor of 
about 2* Thus photoelectron statistics lead to an energy 
measurement uncertainty 


S£ 

where SE is the energy 
I nter est . In f i g ur e 


I 

■— 

^ \SO ' 

loss in MeV in the scintillator of 
28 we show this relationship of 
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In order to interpret the signals measured by the pm 
tubes it IS necessary to have a precise map of I ight 
collection efficiency vs. position in ea'ch of the 
scintillators. Such maps will ultimately be obtained from 
heavy ion calibrations at the BEVALAC at the Lawrence 
Berkeley Laboratory. Preliminary maps were made by 
looking at the amplitude variation of, the photopeaks 
obtained using a col I imated B\ ^-ray source to 
illuminate selected areas of the crystal. Figure shows 
light collection contours obtained in this manner. It is 
emphasized that although this map was obtained in a test 
configuration prior to assembly of the fl ight stacks we 
believe that the essential characteristics of this map are 
similar to those of the crystals in the fl ight stack. It 
can be seen from the figure that over more than "|5^ of +f>e 
crystal area the I ight col lection gradients are <Q.P^/mm, 
and increase to A.»0.5^/mm at the sides transverse to the 
light pipe's viewing direction. 

Since flight data must be corrected using light 
collection maps made during an acceierator calibration, 
and since the cal ibration and the fl ight may occur 
anywhere from several months to more than a yegr apart, it 
IS essential that the response map not change over 
reasonably long periods of time. It is common in this 
type of experiment to use two or more photomultipliers to 
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F_I_GURE 

Relative light collection 
contours measured using a 
y-ray source for sc 

in a test set-up Light 
gradients are also shown 


e f f 1 c 1 ency 
col I I mated 
I nt I I I ator 
col lection 
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view a single scintillator and to sum the resulting 
charges at the input of a charge sensitive amplifier. 
However, as will be discussed below, photomultipliers 
exhibit gain variations with temperature. Furthermore, 
the detailed temperature characteristics of pm tubes vary 
between tubes of identical construction and even vary with 
time for any particular pm tube. Thus the relative 
weighting of the contributions from each tube to the 
overai I response can change in an unknown manner and 
thereby alter the light collection map. In order to avoid 
this problem we have chosen to have only one pm tube 
viewing each scintillator. Gain changes in the tube will 
not affect the light collection of one area of the 
scintillator relative to another but rather will just 
cause an overall normalization error in the pulse heights. 
Techniques by which the normalization can be obtained 
during flight will be discussed later. 

Since the light collection efficiency varies from 
point to point in the crystals, uncertainties in the 
energy loss measurement are produced due to uncertainties 
in the determination of a particle’s trajectory. Light 
collection gradients of combined with a position 
uncertainty of 1 mm yield 

C7FeO> Collection) ^ Q 

5E 
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The pulse height analyzer circuitry introduces 
several sources of uncertainty into our determination of 
the energy deposited in a scintillator. Random electronic 
noise in the linear amplifier stages is responsible for an 
uncertainty equivalent to an energy loss of Q.Qii MeV in 
the low range of the pulse height analyzer and to an 
energy loss of I.l MeV in the high range, so we include a 
relative uncertainty of ^ 



woise) 




0.011 MaV 

1.1 MeV 
Se. 




in our calculation of the overall energy measurement 
uncertainty. The various sources of background which 
stimulate the pm tubes (primarily the •^-ray calibration 
sources) cause the early stages of the pulse height 
analyzers to be continually active at a low signal level. 
Although these pulses are too small to trigger the 
d I scr I m I nators they are capable of causing small base line 
shifts in the amplifiers. For the 0"| scintillator, which 
IS subject to the largest background, we measure 
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The ana I og-to-d I g I ta I conversion introduces uncertainty in 
the energy measurement due to the finite channel width. 



Our logarithmic analyzers have a channel width equal to 
of fli® pulse height being measured, so we obtain 


cn 
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0.00057 


There are a number of systematic effects which will 
introduce significant errors if they cannot be eliminated 
by appropriate calibrations. The gam of the 

photomultipl lers is subject to temperature variations as 
large as -0-8?/^C. However the pm tube gam is not a 
strictly reproducible function of temperature, so a 
knowledge of the pm tube temperature is not sufficient for 
making corrections to the pm tube gain. The scintillation 
efficiency of the Csl(T!) crystals is also subject to 
temperature variations estimated to be (Birks, 

The gam of the pulse height analyzers also 
exhibits a temperature dependence up to 


Three types of in-flight calibrations are planned for 
purposes of correcting for these gam changes. Periodic 
pulser calibrations of the pulse height analyzers at a 
number of amplitudes should allow us to make corrections 
for any electronically prod uced ' var i at i o ns in the signal 
sizes. In order to cal ibrate the photomultipl ler gains a 
calibration source consisting of a small piece of Cs!(TI) 
and a 2 /<C i Cs source has been coupled directly to the 
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face of each pm tube. The logic is designed to 

periodically record non-co i nc i dent events, primarily 
■y-rays. Approximately 3000 ^-ray pulse heights will be 
recorded per hour from each of the eight pm tubes. 

Fitting techniques should allow us to determine changes in 
the pm tube gains with a precision of a few tenths of a 
per cent about once per hour. 

If such a calibration is to be sufficient it is 
essential that any variations of the pm tube temperature 
on time scales shorter than one hour be smooth. To 

produce this condition the pm tubes have been enclosed in 

thermally insulating boxes made of 2" thick slabs of 
polyethylene foam. Two thermistors have been mounted in 
contact with each pm tube to monitor the temperature 
variations seen by the tube. Figure 30 is a plot of the 
temperature measured at one of the pm tubes during a 
balloon flight of 12i. hours duration at float, launched 
from Aberdeen, South Dakota on b June 1311* The variation 
of the temperature was smooth and the maximum temperature 
gradient reached was 0-51®C/hour. This gradient combined 
with a photomultiplier temperature coefficient of 
yields a variation of gain with time of Q.5?/hour. 



FIGUR E 


Temperature variation of pm tubes during 
a balloon flight on b June iSlI* The 
maximum temperature gradient is 


indicated. 







Additional information on the variation of the 
scintillation efficiency and photomultiplier gains can be 
obtained from cosmic ray data itself. The logic assigns 
penetrating events--those which trigger the D8 
scTntillator (primarily relativistic carbon and oxygen)--a 
lower priority than stopping events. These events are 


read out 

only 

i f no stopp 1 ng 

event 

i s 

available 

for 

readout . 

The 

expected 

penetrat i ng 

rate 

(b 

^iO) 

1 s 

approx i mate 1 y 

i500/hour . 

Landau 

f 1 uctuat i ons 1 ead 

to 

an 


uncertainty of in the energy loss of a relativistic 

carbon nucleus in the scintillators. So the mean signal 
produced by such particles should be obtained to 
approximately once per hour. In estimating the 
uncertainty in our energy measurements we assunrie a rms 
gain uncertainty of Q.i^ due to the limited precision of 
our gain calibration resulting from finite counting 
statistics, thus 

0 00 ] 
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It can be seen in figure 28 that at low energy losses 
(^25 MeV) the energy measurement uncertainty is dominated 
by base line shift effects, at energies between 25 MeV 
and ^^i.5 GeV the uncertainty due to photoelectron 
statistics dominates, and at the energies above 1*5 GeV 
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the uncertainty due to light collection variation with 
position dominates. 

The energy measurement uncertainties described above 
apply separately to energy measurements in each of the 
scintillators through which a particle passes. So when we 
combine the energy loss measurements from several 
detectors to form AE we must add the uncertainties in 
these measurements in quadrature to obtain the uncertainty 
in JiE . For those contributions to the energy measurement 

t 

uncertainty which are proportional to the measured energy 
this compounding will reduce the magnitude of the 
contribution relative to' that which would have been 
obtained in a single thick detector. For contributions 
proportional to the resulting uncertainty will be 
unchanged and for those which are independent of energy 
the uncertainty will be increased. 

a /4 

In figure 3i w© show for the isotopes 'Be, 0, Si, 
and Fe the mass resolution as a function of the depth to 
which the particle penetrates in the Csl(TI) stack. These 
results are obtained by evaluating equation 5 using the 
partial derivatives in Table 1 and the uncertainties 
discussed in the preceding paragraphs. A 15^ angle of 
incidence has been assumed in these calculations. In 
addition the contributions to the mass resolution due to 
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FIGURE 21 

Mass resolution calculated for 3e, 0, 


**Si and 

Si>r. 

Fe 

1 nc 1 dent 

at 15* as a 

function 

of 

depth 

of 

penetr at i on 

(R/ secS) 

i nto 

the 

sc 1 nt 1 

1 1 1 ator stack - 


Contributions due to uncertainties in 

measuring AE , E* and L are shown, as 

well as the contribution due to Landau 

fluctuations. The mass resolution 

required to separate adjacent isotopes 

with various relative abundances is 

indicated. The discontinuities in the 
9 

Be curves are due to electronic noise 
near the high level discriminator 
threshold (see figure 28^* 
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Landau fluctuations, thickness uncertainties, and 
uncertainties in the measurements of AE and E* are shown. 
Tick marks on the right hand sides of these plots indicate 
the mass resolution required to obtain an inflection point 
between Gaussian distributions whose centers are separated 
by 1 AMU for various abundance ratios. 
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3 . CALCULATIONS 

Measurements of the isotopic composition of the 
cosmic rays are expected to place significant constraints 
on theories of the synthesis, acceleration and propagation 
of the cosmic rays. Such measurements should make it 
possible to refine estimates of astroph ys 1 ca I parameters 
which have been obtained from elemental abundances and 
should allow determination of additional parameters which 
cannot be obtained from elemental abundances alone. In 
section 3-1 we will briefly review the various classes of 
information which can be obtained from measurements of 
cosmic ray isotopic abundances. More detailed discussions 
of potential applications of isotopic composition 
measurements have been published by a number of authors 
(Meyer, 1335 j Raisbeck et al., 1315b; Shapiro and 
Silberberg, 1315b; Silberberg et al., 131b; Woosley, 
111b). In section 3*2 we discuss the "leaky box model" of 
cosmic ray propagation and the assumptions upon which it 
is based. In section 3*3 we show how in the context of 
this model one can obtain cosmic ray source abundance 
ratios of isotopes of individual elements. We go on to 
evaluate the uncertainties in the ratios caused by 
uncertainties in various measurable parameters and we 
present curves which can be used to directly estimate 
source ratios and their uncertainties in several cases. 
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3.1 Information from Isotopic Composition 

3.1.1 Source Abundances 

Theories of nucleosynthesis predict the relative 
abundances of large numbers of nucl ides which will be 
synthesized under particular astrophys 1 ca I conditions. 
From measurements of these abundances it Is possible to 
derive values for important astrophys 1 ca I parameters. In 
the case of the cosmic rays, however, abundance 
measurements have been largely restricted to the 


separation of 

individual 

elements, due 

to experimental 

difficulties in 

r eso 1 V 1 ng 

isotopes. A 

great deal 

of 

1 n format 1 0 n is 

1 ost 1 f 

only elemental abundances can 

b e 

measured. First 

, when the 

number of observed species 

1 s 

red uced , one's 

ab I 1 1 ty 

to stringently 

constr a 1 n 

the 

n uc 1 eosyn thes 1 s 

theor 1 es 

is d 1 m 1 n 1 shed 

Second , 

the 

influence of 

cer ta in 

astrophys 1 cal 

par ameter S-- 

for 


example, the neutron excess in the source--is strongly 
reflected in the synthesized isotopic abundances but only 
weakly in the elemental abundances. Consequently 
information concerning such parameters is largely obscured 
if only elemental abundances can be measured Woosley 
^lS1bJ has investigated the importance of measuring 


isotopic ratios 

among 

the 

iron peak elements 

for 

determ 1 n 1 ng the 

degree 

0 f 

neutron enrichment of 

the 
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material from which these species were synthesized. 

Furthermore, the elemental abundances observed at 
earth can be obtained from a wide range of isotopic source 
compositions. Therefore, it is necessary to measure the 
local isotopic composition in order to deduce the isotopic 
composition at the sources. In fact, there are cases 
where the isotopic composition must be measured In order 
to obtain even the elemental composition at the sources. 
For example, in a cosmic ray source composed of solar-1 ike 

■k 

material, the calcium will be composed almost entirely of 
^^Ca. However the spallation of ^^Fe would produce 
substantial amounts of secondary ^*Ca, "^Ca, and ^^Ca 

as well as . As a result, the observed calcium would 

be approximately half primary and half secondary in 
origin. On the other hand, if one can observe the isotope 
Ca alone, then one will be dealing with a 
pr imary-to-secondary ratio in excess of 

Additional motivation for the measurement of the 

isotopic composition of the cosmic rays has been provided 

by the suggestion of several authors that the cosmic rays 

are subject to elemental separation effects prior to being 

✓ 

accelerated. Evidence has been presented (Cas.se et al., 
1*3153J Casse and Goret, I311h) which suggests that the 
differences in elemental composition between the cosmic 
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rays and solar system material could be predominantly due 
to such selective acceleration effects. The relative 
abundances of the various isotopes of a particular element 
will not be affected by such charge dependent selection 
mechanisms. It is therefore desirable to use measurements 
of the isotopic composition of individual elements as a 
source of information which will not be altered by 
charge-dependent separation effects. 

3.1.3 Secondary Production and Propagation 

Much of the present knowledge regarding the 
propagation of the cosmic rays through the interstellar 
medium has been derived from the abundances of those 
elements which are believed to be absent at the cosmic ray 
sources. Due to lack of detailed information concerning 
the interaction of the cosmic rays with the 
electromagnetic fields in the propagation medium and of 
the boundary conditions appropriate for describing this 
medium, the propagation of cosmic rays has general ly been 
treated by assuming a distribution of path-lengths through 
which the particles have passed during transport from 
their sources to the point of observation. If the 
path-length distributions described by such a model can be 
characterized by a single parameter (for example, the 
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thickness of a slab or the characteristic length of an 
exponential distribution) then It is possible to evaluate 
this parameter by using available nuclear fragmentation 
cross sections and the observed abundances of a single 
purely secondary species and of its progenitors. Once the 
path-length distribution is known, it becomes possible to 
calculate the secondary contributions to other species and 
then to obtain the primary contribution by subtracting the 
secondaries from the observed abundance. The source 
abundances can then be obtained by applying corrections 
for the nuclear destruction of these particles in passing 
through interstellar matter. 

If, in applying this technique to determine the 
primary component of a particular element, one must use 
only elemental abundances, then it may be necessary to 
utilize an element relatively far from the element of 
interest as the tracer of secondary production. The 
elements below iron which are believed to be entirely of 
secondary origin are the light elements Li, Be, and B and 
the elements Sc, V and Mn in the iron group, plus the 
elements F and Cl. There are no elements which are 
obviously of purely secondary origin near some of the 
elements which one would like to investigate, for example 
Mg. A path-length derived from the light elements will 
reflect the spallation of carbon and oxygen, while a 
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path— length derived from the iron group elements will 
reflect the spallation of Iron. The secondary magnesium, 
however, is produced predominantly by spal lation of SI and 
S. So one is forced to make an additional assumpt i on--for 
example, that the distribution of path-lengths Is the same 
for all species--in order to obtain abundances of many 
elements at the cosmic ray sources. 

in cases where isotopic abundances can be measured it 
IS possible to avoid this difficulty. Many of the even— Z 
elements have one or more isotopes which should, be of 
purely secondary origin. By using a secondary isotope of 
the element for which we want to obtain source abundances 
as a tracer of the secondary production during 
propagation, we obtain a correction which is appropriate 
for this element. 

3.i*3 Radioactive Species 

Various information can be obtained from radioactive 

species in the cosmic rays (Soutoul et al., 1935; 

Raisbeck et al., X9”15a, 1915^). isotopes which ^-decay 

6 7 

with half lives '^10 -10 years can be used as a measure 
of the time since their production. If these isotopes are 
not produced at the cosmic ray sources one obtains a 
measure of the residence time of cosmic rays in the galaxy 
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(as determined by the rate of loss due to nuclear 
destruction and the escape from the galaxy). From 
|S-unstable species which are predominantly of primary 
origin one obtains a measure of the time since synthesis 
of these nuclei in the sources. Unfortunately In the 
region Z ^ 28 nriost ^-decay nuclei are not expected to be 
produced in the cosmic ray sources. Possible exceptions 
are ^®Fe and i (Waddington, 1115). However, these 
species have rather short half lives compared to present 
estimates (Hagen et al., iS11j Garcia-Munoz et al., 
illlb; Muller et al., 1111; Webber et al., 1111) of the 
residence time of cosmic rays in the galaxy of ^10^ years 
(for ^^Fe '?^=3xl0^ yr and for i is estimated to be 
'^2x10^ yr (Casse, 1113) when the absence of orbital 
electrons precludes the possibility of electron capture 
decay). Thus these nuclei should only be of value above 
several GeV/nuc where time di lation will extend the 


effect 1 ve life time 

to 

<«10^ 

years . 




Nuclei which can 

on 1 y 

decay 

by 

means 

of orbital 

electron capture 

can 

be 

used 

as 

probes 

of var 1 ous 


astrophys I ca I conditions. Those electron capture nuclides 


which are synthesized in the source region can be used to 
measure the time between synthesis and acceleration to 
relativistic energies (Casse and Soutoul, 1115^; Soutoul 
et al., 1115; Soutoul et al., HID since only during 
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this time are they expected to have orbital electrons 
available for capture. Electron capture nuclei which are 
entirely of secondary origin will probe the density in the 
propagation medium since after picking up an orbital 
electron from the interstellar gas the nucleus will decay 
only if a time comparable to its electron capture half 
life can elapse before encountering an atom of the gas and 
undergoing an interaction which will re-strip the electron 
from the nucleus. 

In most cases, measurements of elemental abundances 
are of little use in addressing these questions of the 
survival or decay of radioactive isotopes. Other isotopes 
of the same elements generally have significant abundances 
and the difference between complete survival and complete 
decay of the radioactive isotope may only be a few tenths 
of the observed elemental abundance. Uncertainties in the 
abundances of parent nuclei, in nuclear fragmentation 
cross sections and in estimates of the path-length in the 
interstellar medium can easily obscure differences of this 
magn i tude . 



3*2 Cosmic Ray Propagation -- Leaky Box Model 

As mentioned above, various simplified models of the 
propagation of cosmic rays between their production and 
their arrival at earth have been used for unfolding the 
cosmic ray source composition from the effects of 
spallation during propagation. Of these models, the 
•'leaky box model" (introduced by Cowsik et a!., iSb”|) is 
particularly appealing both because it can be obtained 
from a transport equation by invoking a number of simple, 
physical ly reasonable assumptions and' because it has been 
applied with considerable success to observations of 
elemental abundances. In this section we will briefly 
review the leaky box model and the assumptions on which it 
I s based . 

The propagation of cosmic rays can be described by 
the transport equation (Gloeckler and Jokipii, i%S» 
Meneguzzi et al., 1311) 



In this equation n.(r,E,t) is the number per unit volume 
per unit energy per nucleon of particles of species i at 
position r, energy per nucleon E, and time t. The left 
hand side describes local changes in this density due to 
change of position and energy, while the right hand side 
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describes changes of a non-local character (i.e., not 
described by a differential operator) which are due to 
creation (C*, ) and destruction (D* ) of particles. The 
creation term, Cj, is used to describe the production of 
particles of species i in the cosmic ray sources and by 
spallation of heavier species, while the destruction term, 
Dj , describes the destruction of particles by inelastic 
nuclear interactions with the interstellar gas and by 
radioactive decay. This equation assumes that the spatial 
transport of cosmic rays is a diffusive process 
(presumably resulting from the scattering of particles by 
I r regu I ar i t I es in the galactic magnetic field) and that 
the cosmic rays are isotrop.ic. 

Due to our lack of knowledge concerning the spatial 
distribution of cosmic ray sources and of interstellar 
magnetic fields and gas we assume , the diffusion 
coefficient, K , the rate of energy loss, (dE/dt)* , and 
the rates of particle creation, C; , and destruction, Oj , 
are independent of position. Then the spatial dependence 

* Measurements of cosmic ray anisotropies at energies 
of 10“ eV yield values <2x10'"'^ ( see , for example, Elliot et 
al*» 1410>- At lower energies it is only possible to set 
an upper limit on the galactic cosmic ray anisotropy from 
measurements made near earth since interactions with the 
solar wind can produce local anisotropies in excess of the 
galactic value. 
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of nj can only artse from the boundary conditions on this 
equation, which are yet to be specified. In the leaky box 
model it is assumed that the boundary of the confinement 
volume acts as a surface which reflects most of the 
particles incident upon it, but has a small probability of 
1 ett ing a particle escape 'from the volume. 


It is conventional to make the additional assumption 
that the diffusion term in the transport equation can be 
neglected. In this case nj is treated as a constant 
throughout the confinement region and the loss of 
particles at the boundary is accounted for by the 
inclusion of a position independent probability of loss 
per unit time in the term Dj . By making these 
approximations one obtains 
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This approximation can be justified if the cosmic 
rays move freely about the confinement volume, 
encountering its surface many times before finally 
escaping. The typical amount of material traversed by 
cosmic rays before escaping is /^5*5 g/cm^, as deduced 
from measurements of the abundances of purely secondary 
cosmic rays relative to their primary progenitors (for 
example, L i +Be + B/C + N+0 ) . For an Interstellar density of J. 



H atom and 0-i He atoms per cm®, this amount of material 
corresponds to a distance of 80Q kpc (approximately iQO 
times the distance between the sun and the galactic 
center). Since the various regions suggested for the 
confinement of the cosmic rays have scale sizes on the 
order of the galactic radius or smaller, it should be 
possible to neglect spatial gradients if the cosmic rays 
are not strongly scattered while traversing the 
confinement region. We should note that there are 
ind'ications (Cesarsky et al., 1811^ that there may be a 
significant gradient of cosmic ray density in our galaxy, 
if these results are borne out by further observations, it 
will be necessar,y to consider the confinement of cosmic 
rays in a region small compared to the size of the galaxy 
if the leaky box model is to be retained in its present 
f ormu I at I on . 


Next we assume that a steady state has been 

established so that we can set to give 
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The major energy changing process affecting cosmic 


ray nuclei is ionization energy loss due to collision with 
electrons in the interstellar gas. Thus we write 
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where ^ is the density (g/cm^) of the interstel lar gas and 
the other symbols have been defined in section 2 of this 
thesis. Note that although ionization energy loss is 
actually the result of a large number of discrete 
collisions, we can approximate it by a continuous rate of 
energy loss since other relevant length scales are much 
larger than the typical distance between collisions wi+h 
interstellar electrons. For example, a particle with y=2 
(E/M=0|3i MeV/nuc) loses of 'fs energy in collisions 
which transfer less than 150 l^eV to an interstellar 
electron {Rossi, iS52^- The particle's rate of specific 
ionization will be no less than MeV/(g/cm^) (the 
minimum dE/dx for a proton), so on the average it will 
travel no farther than 35 mg/cm^ before losing J.50 KeV . 
This distance is to be compared with nuclear interaction 
lengths which range between 2 g g/cm^ and the escape 
mean free path of '^5-5 g/cm^. Note also that the maximum 
energy which a heavy particle with ^=2 can transfer to a 
free electron is rv3 MeV. Such collisions are rare, but 
energy losses of this size will not significantly alter 
either the ionization rate or the nuclear interaction 
cross sections for particles at the energies being 
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The term, 0- , which accounts for the destruction of 
particles of species i can be written as (Meneguzzi et 

al., mi) 

where n^^ and are the number densities of H and He, 
respectively, in the interstellar medium, and <T^i are 

the cross sections for breakup of species i in a collision 
with an interstellar H or He nucleus, respectively, 
is the mean time for escape from the confinement volume, 
and 'T] is the mean time for radioactive decay in the frame 
in which the particle is at rest for stable 

species). The three terms on the right hand side of 
equation "1 represent loss by means of nuclear destruction, 
escape from the confinement region and radioactive decay, 
reading from left to right. 

The creation term, C,* , is written 

J 

where is the rate of production (per unit volume, per 

unit energy per nucleon) of cosmic rays of species i by 
the cosmic ray sources and C^j,(E,E') and j ( E , E ' ) are 

the cross sections for species j at at energy per nucleon 
E* to produce species i at energy per nucleon E in a 
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collision with interstellar H or He, respectively. There 
IS evidence (Greiner et al., 1415^ that at the energies of 
interest here (SlOO MeV/nuc) spallation products are 
predominantly produced with velocity (or energy per 
nucleon) nearly equal to that of the incident heavy 
nucleus. Therefore we obtain 


It IS convenient to replace the cosmic ray density n* 
by the omnidirectional flux, jS/. and to define 

q j SQj / ( n^M^ + n^g ) . We obtain 



A 
A esc 




— Of) 







In this equation we have written have 
assumed that we are not dealing with any radioactive 
species, so we can neglect destruction by radioactive 
decay. Equation 8 is the basic equation relating observed 
cosmic ray fluxes ()2^{ ) to source abundances (q,) in the 
leaky box model. 





To simplify notation 


we 


will make the following 


definitions; 


^ 'Ah * (n„yn„) 




gpil'*' (^Me/llM^OTi! 
Mh + (h„,/n„)M.e 


Then we obta i n 


0 = 



3.3 Secondary Tracer Formalism 

We will now specialize to the case of the various 
isotopes of a single element. Thus we will drop the 
subscript on the nuclear charge, Z, in equation 3- Note 
that in equation 3 The j2f. are the cosmic ray fluxes 
measured in local interstellar space and the "secondary 

A 

source term", S. , depends only on these fluxes and on 
nuclear cross sections. Thus the only quantities in 
equation 3 which we cannot directly measure are the source 
abundances, q^ , and the "escape mean free path", Aes<, • If 
we know that a particular species is not produced in the 
cosmic ray sources (qj=Q) then we can obtain A e** 1 n terms 
of measurable quantities (Reames, 131H)t 



A 
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hi M 



If we substitute this value of ^esc obtained from species 
I into equation written for species j, we obtain 



This equation can be simp) ified by approximating 

and s(E) in the vicinity of the energy of interest, E®, by 

power laws: 



Then the ionization energy loss terms reduce to 



We can write the analogous expression for the source 
abundance, q., of a third species, and the ratio of source 

K 
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abundances can be written in the form: 



The right hand side of equation 12 depends only on 
measurable quantities (local interstellar cosmrc ray 
fluxes, spectral indicies and nuclear cross sections). 


Below we will discuss the degree of precision to which 
these quantities must be known in order to obtain useful 
estimates of the ratio of source abundances q./q, . 

i] 

In evaluating equation 12 we employ f r agmentat i on 
cross sections obtained from the sem i -emp i r i ca 1 formulas 
of Silberberg and Tsao (1443a) with modifications given by 
Silberberg et al. Total inelastic cross sections 
are obtained using the formulas (Tsao and Silberberg, 
1415» Cheshire et al., 1414)5 


g-,= St mb (l-atTMr*') 


Notice that in evaluating equation J.2 we require 
ratios of the fluxes of parent species (used in 
calculating , S. and to the fluxes of the isotopes 
under consideration. However, at present, only elemental 
abundance ratios are available. For this reason we shall 
assume that the three isotopes of interest constitute the 
entire local abundance of the element being considered. 
This IS an adequate approximation in all the cases which 
we shall examine if the source abundances are not vastly 
different from the solar system abundances. This 


assumpt i on 

will not be a 1 

imitation 

when 

1 sotop 1 c 

abundance 

observations of the elements 

wh i ch 

we 

are 

considering 

become aval 1 ab 1 e , 

since it 

will 

then 

be 

poss i b I e 

to appropr 1 ate 1 y 

renorma 1 i ze 

the 

e 1 ementa 1 

abundances . 






Note 

that the secondary 

tr acer 

technique 

for 


obtaining source abundance ratios can only be applied in 
the case of elements having an isotope which is believed 
to be absent at the cosmic ray sources. Among such 
elements there are several which are of considerable 
importance for the theory of nucleosynthesis, for example 
S, Ar and Ca (Casse and Meyer, i'^Ha). 
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3.3*1 Analysis of Uncertainties — Motivation 


Before beginning a detailed treatment of the effect 
of uncertainties in the measured quantities in equation 
on the deduced source ratio we wish to motivate this 
treatment by means of a sem i -quant i tat i ve discussion of a 
simplified case. For this purpose we shall make the 
following assumptions; 1) ionization energy loss is 
negligible (s(E) = f))> 2^ the total inelastic cross sections 
for the species of interest are equal ^ A;jp, = » 3^ 

isotope k is significantly more abundant in the local 
cosmic rays than species j or i ), 4) +he 

observed flux of species j is largely secondary in origin 
while that of species k is largely primary in origin, and 
5) cross sections for species i, j and k to make one 
another are negligible. Using assumptions (1) and (3) we 
can reduce equation 13 to 


% 



Assumption (4) allows us to write 


k = 3L 

Ik A 



This expression illustrates the difficulties involved in 
obtaining precise values of cosmic ray source ratios. The 
assumption that the observed flux of species j is largely 



i5i 


of secondary origin implies that the difference of the two 
terms on the right hand side of this equation is 
significantly less than either of these two terms 
individually. Therefore a small relative error i rt either 
of these two terms will translate into a relative error in 
q^/qj^ which is larger by a factor of approximately 
02<-/>2^h)/(q/qk)- 


In order to make this discussion more precise, we can 
calculate the uncertainty in the calculated source ratio 
due to uncertainties in the measurement of the local 
fluxes >2^. and and in the determination of the 

A A 

secondary correction factor, R3(S./Sj). To simplify the 
notation we will write r^^ “ ’ 
Our assumptions (3) and (5) allow us to treat r.^^ , r,|^ and 
R as statistically independent. We obtain 





If we now use assumptions (3) and (H) and also assume that 
the uncertainties in the observed fluxes are entirely 
statistical and given by k ' ^ ^ 

wr i te 



where N is the observed number of events of species k (and 
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approximately equal to the combined total number of events 
of species i, j and k). This expression explicitly shows 
the magnification of both the statistical errors in the 
measured fluxes and the errors in the calculated quantity 
R by a factor of r.,/o(. which is likely to be >5 for 

J H J 

most of the cases examined below. 


3-3*2 Source Ratio Uncertainties 

We now wish to consider in detail the effects of 
uncertainties in the measurable quantities — local fluxes, 
spectra and nuclear cross sections--on the uncertainty in 
the source abundances that we derive using the secondary 
tracer technique. By so doing it is possible to identify 
those areas in which more precise measurements can improve 
our ability to derive cosmic ray source abundances. 

For purposes of this discussion we restrict our 
attention to a particular set of possible source 
abundances. In particular, we use the elemental source 
ratios obtained by Silberberg et al. (13%) from 
available cosmic ray elemental abundance measurements near 
earth and from the assumption that the propagation is 
described by an exponential distribution of path- I engths . 
We follow their procedure (Tsao et al., 1S13» Shapiro et 

1315a) of distributing the source abundance of each 
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element among its stable isotopes in proportion to the 
solar system abundances compiled by Cameron We 
then set the source abundance of the istotopes. which we 
are using as tracers to zero. Table 1 lists the resulting 
source abundances. For the tracer isotopes we indicate in 
parentheses the source abundances which we had before 
setting these abundances to zero. 

The source abundances are then propagated through an 
exponential path-length distribution with g/cm^ 
in a medium composed of H and He in a ratio He/H=n*i by 
number of atoms. This calculation was performed for an 
Interstellar energy of HGD MeV/nuc and ionization energy 
losses were neglected. Table 1 also shows local 
abundances obtained using this procedure. 

In table 8 we summarize the effects of variou? 
uncertainties on the source ratio predictions obtained by 
applying the secondary tracer technique to isotopes of the 
elements S, Ar, and Ca. The local abundances 
listed in table ”| are employed as a basis for calculating 
source abundances using equation 12? Since the local 
abundances in table were calculated neglecting 
ionization energy loss, when equation 12 is applied using 
these local abundances it should only reproduce the source 
abundances (also shown in the table) if we again assume 
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Table 2 

Source and Local Abundances* 


1 sotope 

Source 

Loca 1 

Isotope, 

Source 

Loca 1 

Nr-G’'r 

Tm" 

.00412 

Ar-W " 

O' ■ 

.a‘T57 

Ni-b2 

.03 

.0142 

Ar-38 

-13 

• bbl 

Ni-bl 

• 01 

.00548 

Ar-31 

0 

-251 

Ni-bO 

.21 

.101 

Ar-3b 

• b1 

.148 

Ni-54 

0 

.00501 

ci-31 

.20 

.300 

Ni-58 

•54 

.2bl 

ci-35 

.08 

• 4b1 

Ni-5b 

0 

.0D13b 

s-3b 

.0004 

-050b 

Co-54 

0 

• 00352 

s-34 

.13 

-514 

Co-51 

0 

.0208 

s-33 

0 (.02) 

.344 

Fe-58 

.01 

.0311 

s-32 

2.85 

2.15 

Fe-51 

.48 

.242 

p-31 

• 20 

-b20 

Fe-5b 

20.11 

4.8b 

Si-30 

• b5 

.403 

Fe-55 

0 

.511 

Si-24 

• 44 

1-08 

Fe-54 

1.28 

.405 

Si-28 

14-3b 

13-0 

Mn-55 

.10 

.331 

AI-21 

2.30 

2-84 

Mn~54 

0 

.331 

A 1 -2b 

0 

.528 

Mn~53 

0 

.521 

Mg-2b 

2.b8 

2-38 

Cr-54 

.01 

.DbOb 

Mg-25 

2.43 

2.82 

Cr-53 

.04 

.234 

Mg-24 

18-84 

14.1 

Cr-58 

.33 

1-04 

Na-23 

• 40 

2-11 

cr-5i 

0 

.544 

Ne-22 

l-b3 

2.81 

Cr-50 

.02 

.341 

Ne-21 

0 (-04) 

1-51 

v-51 

0 

• lib 

Ne-20 

13-33 

11.4 

V'50 

0 

.281 

F-14 

0 

2.15 

v-44 

0 

.5bb 

0-18 

.23 

1.48 

Ti-50 

.005 

• 03b4 

0-11 

0 (-04) 

1.84 

Ti-44 

.00b 

• 104 

0-lb 

110.1 

84. Q1 

Ti-48 

.01 

.114 

n-15 

-03 

11.04 

Ti-41 

.001 

■ bbl 

N-14 

1.41 

13-84 

Ti-4b 

.008 

.588 

c-13 

Ml 

b.QO 

Ti-H4 

0 

-0308 

c-12 

48-84 

44-00 

Sc-45 

0 

.484 

B-11 

0 

11-b 

Ca-48 

.004 

.00225 

B-10 

0 

8-13 

Ca-Hb 

.00008 

•0011b 

Be-10 

0 

l.b4 

Ca-44 

• 05 

• 41b 

Be-4 

0 

2-80 

Ca-43 

0 (.003) 

■ 343 

Be-1 

0 

5-04 

Ca-H2 

-02 

• 4b2 

Li-1 

0 

b-84 

Ca-4i 

0 

.104 

L 1 -b 

0 

8-31 

Ca-40 

2.33 

1.34 




K-Hi 

• O01 

.240 




K-40 

0 

.242 




k-34 

.04 

• 44b 





* Normalized so that carbon=lQ0. 
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that s(E)=Q. In line 1 of table 8 we show the source 
abundance ratios from table 1 ( *^Ne/*^Ne, 

^^Ar/ ^Sr, *^^Ca/*^^Ca and *^^Ca/^‘’Ca) for the cases 
which we are examining. Line 2 shows the values of these 
ratios which are obtained by applying equation [Z with 
s(E)=0 to the local abundances listed in table 1. The 
agreement between lines 1 and 2 's not perfect, the 
residual differences of ^2% being artifacts of the two 
numerical calculations. These differences are small 
compared to the important effects which we wish to 
exam i ne . 

Having investigated our ability to reproduce sourjce 
abundance ratios using equation 12 we now wish to include 
the effects of ionization energy loss, since only when we 
do so will we be able to investigate the effects of 
uncertainties In the spectral shapes of the species of 
interest. In line 3 we show the source ratios resulting 
from the re-eva 1 uat i on of equation with the inclusion 
of the correct values of s(E) (at IQQ MeV/nuc, 
s(E)/E— .00b cm^/g) but still using the local abundances 
shown In table The differences in the values of the 

source ratios calculated with and without the effects of 
ionization energy loss range up to 18? in the cases 
examined. Clearly the approximation that s(E)==0 should be 
avoided even at energies as large as 300 MeV/nuc. In the 



Tab le \ 

Sources of Error In Calculated Source Ratios 



Primary Isotopes 
Tracer Isotope 
Descr t pt I on 
J. ) Cameron qj/qj, 

equation 12 with s(E)=0 
3) equation 12 with 
correct s(E) 

H) parent Isotopes varied 


'‘O, '*0 

•»o 


OHxlT* 


1 3Hxl{T= 


rms deviation resulting from 


5> variation of parent 
I sotopes 

b) variation of parent 
el ements 

I) 20? error in fragmentation 
cross section 

norma! i zatlon 

8a) 35? uncorrelated errors 
in fragmentation 
cross sections 

8b) linear approximation 

8) 35? errors In 

fragmentation cross 
sections — parti y 
corre I atad 

10) 10? errors in total 
Inelastic cross 
section normalization 

II) 10? uncorrelated errors 
In total inelastic 
cross sections 

12) errors of 0-2 In 
spectral indicias 

13) error due to statistics 

in measuring and 

0 /j 2 t using 1000 events 


0 2 >«ifr^ 
0 2 xio’^ 
0 1 ><1C'* 

H 0 xiCT* 



0.1 xiO'"’ 

0 1 xiO'* 

0 H xlO'^ 

1 1 xlO'* 


‘*"Ne, 

^'Ne 

j: 

‘s, ”s 

^*Ar . ”Ar 
”Ar 

"Ar. ”Ar 
»‘’Ar 

‘’"Ca, **Ca 


’Ca. ’’Ca 
■*^Ca 

1 223xlo;: 
1 232x10 

4 

bGxlT"* 

1 

440x10. 

i 

040x10 

8 bxin'* 

2 

15x10'* 

4 

SbxiO'*’- 

1 

434x10’' 

i 

030x10" 

8 5x10'* 

2 

14x10'* 

L.20bxlf)-' 

4 

3bxl0"^ 

i 

424x10“' 

2 

OlbxiO” 

10 1x10’* 

i 

01x10'* 

1 185xiT’ 

4 

13x10"^ 

0. 

332x10" 

4 

500x10" 

1 8x10'^ 

4 

08x10’* 

0 OOSxiO*' 

0 

54x10'* 

0 

40 xlO"' 

0 

55 xlo" 

0 1x10'^ 

1 

1 xlO'*' 

0 OHGxiOr' 

0 

51x10'* 

0 

50 XlO” 

0 

15 xin" 

4 bxlO"’ 

n 

82x10"* 

0 01 xiQ“ 

0 

05x10"* 

0 

025x10"' 

0 

022x10" 

0 4x10"* 

0 

04x10’* 

n 31 xio" 

4 

0 xl0“* 

4 

b xio" 

1 

4 xlo" 

b4 xlO'* 

[ 

) bxin’* 

0 33 xiQ-; 
0 10 xiO" 

4 

i xlO'* 

2 

1 xir' 

? 

3 xlo" 
18 xlO" 

b3 xiO'* 

b 

1 xlO"* 

1 

1 xlO"* 

i 

2 xlQ” 

0 

4 3xlT* 

0 

80x10'* 

% 

0 nnixio" 

0 

04x10'* 

0 

022x10' 

0 

OlbxlO" 

0 Ixlif* 

(1 

03xiir* 

0.008x10"* 

0 

05x10"'' 

0 

021x10" 

0 

OSbxlO" 

0 5x10'* 

0 

04x1 o'* 

0 055x10"' 

0 

83x10'* 

0 

53 xiO'’ 

0 

41 xlo" 

13 xiT* 

0 

SOxiO'* 

0 28 xl0"‘ 

3 

b xlO'* 

1 

4 xlO" 

1 

0 xiO'' 

31 xiO'^ 

3 

1 xin’* 
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analysis which follows of other sources of error in 

calculating the source abundance ratios we have included 
the effects of ionization energy loss. Therefore the 

resulting ratios, q /q., should be compared with the 

J * 

values shown in line 3 of table 8 rather than with those 
in lines J. and 2* 

In calculating the effects of 'the various sources of 
error we have used a Monte Carlo technique in which the 
parameter being investigated is varied with an appropriate 
distribution and the remaining parameters are held fixed 
at their nominal values. We generate 300 such cases and 
calculate the rms deviation of the values of q /q, 

obtained from equation 13- 


3*3*2*i Abundances of Parent Nuclei 


The secondary source term. 



Is of the form 








A 

where the summation is overall species which 
species i by col I ision with interstel lar 


can produce 
hydrogen or 


helium. Since it is experimentally more difficult to 
resolve the isotopes of elements of higher atomic number 
than those of lower atomic number (see, for example. 
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figure 3i)> < is important to address the question of 
whether detailed knowledge of the isotopic compos f t i on _ of 
parent nuclei is required in evaluating S; , or whether 
the elemental abundances which are presently available 
(Senegas et al., iS15? Ormes et al , 1315; S i I berberg et 
al'» iS1b; Caldwell and Meyer, 1311; Gar c i-a-Munoz et 
al** ISlIa) are adequate. in order to answer this 
question we have assumed that the local abundances of the 
three isotopes of interest are known and have used values 
of these local abundances from table 1- However we also 
assume that only elemental abundances of the heavier 
(parent) species are known. In order to evaluate equation 
12 we must make some assumption about the isotopic make-up 
of the parent elements. We have chosen to assume that 
there is no a priori preference for one isotope of a 
particular parent element over another. Using this 
assumption ws investigate the importance of knowing the 
isotopic composition of the parent elements. In 
evaluating equation 12 we have used nominal values of 
cross sections, spectral indicies and elemental abundances 
while randomly varying the isotopic composition of each 
parent element. The isotopes that are stable in the 
cosmic rays are randomly assigned fractions of the total 
elemental abundance with equal weights. Having assigned 
abundances, we calculate the S, and then evaluate q./q. 



for the flux ratios and obtained from table 1- 
This procedure was repeated 300 times to obtain a mean 
value of q./q. and a standard deviation due to isotopic 
variations in parent nuclei. The resulting means are 
shown in line H of table while the standard deviations 
are shown in line 5 • 

It can be seen that the mean value of (line H 
of table 8^ obtained while varying the isotopic 
composition of the parent elements is in some cases 
significantly different from the value in line 3 oi the 
table, this difference sometimes being several times 
I a r ger than the rms deviation (line 5^ resulting from this 
variation. This is due to the fact that isotopes which 
are dominant in the solar abundances are those having the 
smallest neutron excess, at least up through Z=3P 
(Cameron, iS13) and that spallation reactions which form 
products which have neutron excesses close to that of the 
parent nucleus are more probable than reactions which 
significantly alter the neutron excess. The largest such 
error occurs in the case of the isotopes of Ar . In this 
case the secondary production is largely due to spallation 


of calcium. 

The 

ratios 

of 

cross sections for 

produc 1 ng 

the isotopes 

"Sr, 

^V, 

"*Ar 

and r in a 

p + 

reaction are 

' i 

: i.H : 

i-H 

: 0 while in a 

p + 

react i on they 

are 1 

: 0.3 

s 5 

: 3 ( S I 1 berberg 

and Tsao, 
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Because of this wide variation a knowledge of the 
isotopic abundance of Ca in the cosmic rays is essential 
for interpreting the observed abundances of the Ar 
i sotopes . 

In most cases the dominant contributions to the 
secondary production of an element come from other 
elements with atomic numbers close to that of the element 
of interest. For example, the secondary 0 is largely 
produced from Ne and Mg, the secondary Ne from Mg and Si, 
the secondary S from Ar and Ca, the secondary Ar from Ca 
and Fe, and the secondary Ca from Fe. It should be 
possible for an instrument capable of measuring the 
isotopic composition of the element of interest to also 
obtain measurements of the isotopic composition of these 
dominant parents. The use of such measurements in 

A 

calculating S? should significantly reduce the uncertainty 
(both the rms uncertainty and the systematic shift) i,n the 
secondary corrections required. 

It should be noted that uncertainties in q./q. due to 
uncertainty in isotopic composition of the parent elements 
can be significantly reduced without having to resolve 
individual isotopes of these elements. Mass distributions 
with mass resolutions as large as i AMU should be of 
considerable value, since spallation cross sections tend 
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to vary smoothly with particle mass. In some cases such 
distributions are already available (Fisher et al., iS1b)- 
In the case of Ar where the dominant parent, Ca, has many 
stable Isotopes, it should be particularly useful to 
eliminate some of these isotopes as major contributors to 
the production of secondary Ar, 


In addition to the present lack of information 
concerning the isotopic abundances of the heavier elements 
there are uncertainties in the elemental 

A 

abundances which will. produce uncertainties in the S. and 

hence in q,/cj. . We have used the tabulation by Silberberg 
j * 

et al. (iSlb) of local elemental abundances and their 
uncertainties, based on the measurements of a number of 
investigators at rigidities '^4 GV (energies 
Sic!15 MeV/nuc). Local abundances of the various elements 
were allowed to vary independently with Gaussian 
distributions with standard deviations equal to the 
uncertainties given by Silberberg et al. (ISTb)* These 
uncertainties range between 2% for the most abunda.nt 
elements (such as oxygen) to 50^ elements with small 

abundances (such as scandium). The, isotopic compositions 
of all elements were simultaneously varied in the manner 
described above. Since the isotopic and elemental 


compositions were varied Independently, we obtained the 


uncertainty in q /q. due to elemental composition ' 

j k 
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uncertainties alone by quadratic subtraction of the 
uncertainty from varying only the isotopic composition 
from that obtained while varying both elemental and 
isotopic compositions. In all the cases considered it was 
possible to obtain a statistically significant result by 
using 200 events. The uncertainty obtained in this manner 
IS shown in line b table 0* 

In comparing lines 5 and b o'f table 8 ■'f' can be seen 
that the uncertainty in estimating q,/q due to present 

J K 

uncertainties in elemental abundances is comparable to the 
uncertainty due to our total lack of information 

concerning the isotopic make-up of the various parent 
elements. Thus a detailed knowledge of the isotopic 
composition of these parents will not in most cases 
significantly improve the accuracy of calculations of 
source abundances until the elemental composition can be 
measured with increased precision. However, presently 
operating particle spectrometers have sufficient elemental 
resolution so that increased exposure above the earth’s 
atmosphere would make it possible to reduce the elemental 
composition uncertainties to a point where some knowledge 
of the isotopic composition of parent species will be 
useful . 
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3*3*E*2 Fragmentation Cross Sections 

In order to investigate the effects of uncertainties 
in the fragmentation cross sections, and 0^^, , on the 

calculated source abundances we have varied these cross 
sections in several ways. First, we have assumed that the 
uncertainties in these cross sections are entirely 

systematic, and in particular that there are no 

uncertainties in the ratios of the cross sections to one 
another but that there is simply an overall normalization 
error. One can o.btain estimates of the magnitude of 
possible normalization errors by examining ratios of 
values obtained for the same cross sections from different 
experiments or calculations. Lindstrom et al. (lS15a) 
found an average ratio of i.2? between their measurements 
of cross sections for the production of lighter isotopes 
by the spallation of and ^^0 and the value obtained 

from the sem i -emp i r i ca 1 formulas of Silberberg and Tsao 
(I313a)- In I I ne ~\ of table g we show the uncertainty in 

q /q. obtained due to a 20^ uncertainty in the absolute 

i 

normalization of the fragmentation cross sections. Even 
for a normalization uncertainty of this size, the 
uncertainty produced in q./g is comparable to 

4 k 

uncertainties due to various other effects and at present 
IS not a limitation on our ability to obtain cosmic ray 


source abundances. 
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Second, we have assumed that the uncertainties in the 
various spaiiation cross sections are entirely 

r 

uncorrelated. We independently make a random selection of 
each of these cross sections from a Gaussian distribution 
with mean equal to the cross section calculated using the 
formulas of Silberberg and Tsao and a selected fractional 
standard deviation. Lindstrom et al. obtain a 

standard deviation of 31^ < ri the distribution of the 

ratios of their carbon and oxygen spallation cross 
sections to those calculated using the formulas of 
Silberberg and Tsao. Also, comparisons by Silberberg and 
Tsao 11 * 3133 ) of their values to the measurements on which 
the sem i -emp I r i ca I formulas are based indicate 
uncertainties of ^30“M0/5. The resulting uncertainties in 
the calculated values of q,/q, are shown in line ga of 

J IC 

table 8 fof" 35)^ standard deviations in the fragmentation 
cross sections. 


In most of the cases being considered the 
uncertainties in uncorre I ated cross section 
errors is at least as large as the value of q^/qj^ itself. 
In section 3-3>H we will ex;amine the extent to which these 
errors must be reduced in order to derive meaningful 
parameters character 1 z 1 ng the cosmic ray sources. In that 
calculation we will assume that the source ratio errors 


I n 


scale 


proportion to the uncertainty 


in the 
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fragmentation cross sections. In order to establish the 

validity of such scaling of the errors in q./q. in 

j h 

proportion to the fractional uncorreiated un certa i n i tes in 

the fragmentation cross sections, ’we have recalculated the 

values in line gs of table 8 by assuming cross section 

errors of 3*5^ and then multiplying the resulting error in 

q^/qj^ by a factor of IQ. These errors obtained from this 

"linear approximation" are shown in line gb of table 8- A 

comparison of lines ga and gb shows that linear scaling of 

the errors in for cross section errors less than 35/^ 

is reasonable except in the case of argon. The failure of 

this scaling in the case of argon is due to the fact that 

3 & 

even the most abundant isotope, Ar , has a large 
secondary component, so that the denominator in equation 
13 IS relatively small. The variation of our estimate of 
the secondary correction due to large Gaussian 

errors in the fragmentation cross sections results in a 
non-Gaussian distribution of q^/q^^ values. The long tail 
on this distribution grossly affects the calculation of 
the rms spread of the distribution. 

Third, we have qssumed that ratios of cross sections 
for producing different isotopes from a single parent 
species can be precisely determined but that the absolute 
magnitude of these cross sections is uncertain by 35/^- 
Such a situation should result when fragmentation cross 
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sections for producing various isotopes are measured in a 
single experiment such as those performed by Heckman and 
CO I I aborators (Lindstrom et al., in such 
experiments details of the experimental set-up such as 
target thickness and beam flux will affect the absolute 
cross sections obtained but should only weakly affect 
cross section ratios. In line 9 oi" table 8 we show the 
uncertainty in q /q, produced by a 35/^ error of this type. 
We have also calculated these uncertainties using a linear 
extrapolation from those obtained assuming a 3-5^ cross 
section error. Good agreement between these two 
calculations” is obtained in each of the cases being 
cons I dered . 

A comparison of values in I i ne 8 with those in I ine 8 
shows that an improvement by a factor of between 3 and 10 
in the effect of fragmentation cross section errors on the 
uncertainty in calculating source abundance ratios will 
result if the cross sections from a single parent can be 
obtained with high relative precision. 

In evaluating the effect of cross section 
uncertainties on source ratio calculations we have used 
cross section errors based on comparisons between the 
sem I -emp I r I ca I cross sections of Silberberg and Tsao and 


various measured cross sections 


These 


errors 


are 
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appropriate for the present calculation since the lack of 

measurements of most of the required cross sections forces 

us to utilize the sem i -emp i r i ca I results. However, it 

should be noted that frequently the uncertainties in 

measured fragmentation cross sections are somewhat less 
% 

than the 35? assumed here (Silberberg and Tsao, 1513b; 
Lindstrom et al., i515a? Perron, 151b)* Thus measurement 
of cross sections for those reactions which dominate the 
secondary production of species which we are considering 
should significantly improve the accuracy of the 
calculations. Note also that the secondary production of 
many cosmic ray species may be dominated by contributions 
from a few nearby parents. In this case substantial 


reduction of the uncertainty in q./qj^ may be achieved by 

J 

improved measurements of a few key cross sections. In the 


present study we have not attempted to identify these key 


cross sections. 


3*3*3*3 Total Inelastic Cross Sections 

We have evaluated the contributions to the 

uncertainty in q./q|, which result from uncertainties in 

J 

the total inelastic cross sections, OT. and OZt . As was 
done in the case of fragmentation cross sections, we have 
examined two extreme cases--cross sections subject to a 
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normalization error but no relative error and cross 
sections subject to uncorred ated errors. Comparisons of 
total Inelastic cross sections measured by fragmenting 
beams- of C, 0 and Ar on a variety of targets (Cheshire et 


al,, 1S1H; Li 

ndstrom et 

a I . , 

iS15b) 

w I th 

simple 

geometrical models 

shows that 

such 

models 

yield 

cross 

sections accurate 

to 

We 

have used 

m 

as an 


estimate of both the normalization error and the 

uncorrelated errors in the total inelastic cross sections. 

/ 

The uncertainties in q, /q, resulting from these errors are 

J K 

shown in lines J.Q and H, respectively, of table 8- It 
can be seen that these errors make relatively small 
contributions to the overall uncertainty in q./qj^. 

3-3-2-H spectral Indicies 

Uncertainties in the spectral indicies of the 
observed fluxes of isotopes of the elements of interest 
contribute to the uncertainty in the derived source ratio 
since ionization energy loss continually reduces the 
energy of individual cosmic ray particles resulting in an 
effective sink (or source) of particles at any particular 
energy. The magnitude of this effect depends on the slope 
of the equilibrium spectra. We have assumed that at the 
energies of interest the local interstellar spectra of the 
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species which we are considering are proportional to the 
proton and alpha particle spectra obtained by Garrard 
from near-earth measurements and calculations of 
the effects of solar modulation. Garrard found a spectrum 
of the form 

With m=93i MeV/nuc. He obtains y<^=0 * 35 i 0 • 3 • From this 
spectral shape we obtain the spectral index at 300 MeV/nuc 
as 

V e 1.33 ±0./?, 


It should be recognized that this value and its 
uncertainty can only be treated as an indication of those 
which should be obtainable from near-earth spectra of the 
species of interest if these spectra can be measured with 
greater statistical accuracy than they have been measured 
at this time. In order to reduce the uncertainties' in the 
spectral indicies it will be necessary to refine our 
understanding of the solar modulation process or to make 
abundance measurements outside of the solar cavity. 

In line 12 of table 8 we show the uncertainty in 
q /q^ which results from independent uncertainties of Q.2 
in the spectral indicies of both of the primary isotopes 
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and of the secondary tracer isotope. For the cases which 
we are considering this source of uncertainty will become 
important when the uncertainty due to uncorrelated 
fragmentation cross section errors can be significantly 
reduced. It should be noted that since the spectral 
indicies appear in equation multiplied by s{E)/E, the 
effects of uncertainties in these indicies will decrease 
with increasing energy. 

3*3*2*5 Isotopic Abundances 

Finally, we show in line 13 of table 8 
uncertainty in q;/qt due to uncertainties in measuring the 
isotopic ratios and . We assume in this and in 

all subsequent sections that the uncertainties in these 
ratios are entirely statistical (i.e>,, that measurements 
can be made with good mass resolution). We have 

calculated the uncertainties in q./q^ for the case of 100(3 

J " 

events observed among the three species (?, j, and k). It 
can be seen from table H that for the elements considered 
here the HEIST instrument will accumulate between 100 and 
5000 events per day over its entire energy range. The 

uncertainties shown in table 8 due to uncertainties in the 
measured isotopic ratios are proportional to 1/<\[n*, where N 
IS the number of events observed, so this uncertainty can 
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easily be scaled to the number of events actually observed 
I n an exper i ment . 

The uncertainties in measuring isotopic abundance 
ratios are generally the most significant contributions to 
the uncertainty in the calculated source ratio along with 
the uncertainty due to uncorrelated errors in 
fragmentation cross sections (assuming iOOO events are 
observed). In order to reduce this uncertainty it will be 
necessary to obtain increased exposure factors 
(geometrical factor times duration of exposure) and to 
obtain sufficient mass resolution so that systematic 
errors in the measurement will not limit one’s ability to 
reduce the uncertainty in the measured ratios to the 
stat i st I ca I limits. 


3-3»2*b Summary of Uncertainties 

At present, the precision attainable in deriving 
isotopic source ratios is limited primarily by 
uncertainties In fragmentation cross sections and by the 
1 imited counting statistics obtained by typical cosmic ray 
experiments. The first of these uncertainties can be 
improved by identifying those cross sections which 
dominate the production of the element of interest and 
then making high precision measurements of the ratios of 
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these cross sections. In fact, significant progress has 
been made both at Berkeley (Lindstrom et al., ISlSa) and 
at Orsay (Perron, lS1b; Raisbeck and Yiou, illla; 
Raisbeck et al., illlb) in the past few years and it 
should soon be possible to utilize measured cross sections 
rather than those obtained from the sem i -emp i r i ca ! 
formulas in many cases. The uncertainty due to limited 
counting statistics will be reduced when it becomes 
possible to fly experiments with large geometrical factors 
(SiOO cm sr) for extended periods on spacecraft outside 
of the garth's atmosphere. 

The systematic uncertainty which is present in some 
cases due to uncertainty in the isotopic make-up of parent 
elements can be reduced, as mentioned above, by using low 
resolution mass spectra to establish limits on the 
distribution of parent isotopes. Clearly, even greater 
improvement will result if isotopes of the Important 
parent species can be individually resolved. 


Of POOIt 
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3-3*3 Curves for Interpreting Local Abundances 


We 1 

now wish to extend 

the treatment 

of source 

ratio 

uncerta inties 

in section 

3.3.3 

to a 

range 

of 

focal 

1 sotop i c 

abundances other 

than 

those 

wh i ch 

wou Id be 

obtained 

from 

a so 1 ar- 1 1 ke 

sour ce . 

By so 

do 1 ng 

It 

will be 

poss i b I e 

to 

character i ze 

the 

dependence 

of 

sour ce 


abundance uncertainties on the magnitudes of the abundance 
ratios observed in local interstei lar space. These 
results used in conjunction with m'odels which predict the 
dependence of the source abundances on various 
astrophys I ca i parameters make it possible to evaluate the 
significance of conclusions concerning these parameters 
which are derived from cosmic ray observations. 
Conversely, these results can be utilized for the purpose 
of designing experiments capable of distinguishing between 
alternative astrophys i ca I models at a predetermined 
significance I eve 1 . 


We separately 

consider 

the 

following triplets 

of 

isotopes: ^^ 0 , 

'*'0, '’O; 


^=^Ne, ^'Ne; , 


“^^Ca, *^^a. 

and 


^Ca, ^^Ca. In 

each 


case the last isotope listed is assumed to be absent at 
the cosmic ray source. The isotopic abundances of these 
species can yield information concerning both the 
nucleosynthesis processes responsible for the production 



of the cosmic rays and the state of evolution of the 
material from which the cosmic rays are derived. In 
section 3'3*H w® will discuss how local abundances are 
affected by one of the important source parameters. 


For 

each case we 

have 

ob ta 1 n ed 

three 

f am 1 1 1 e s 

of 

curves . 

First, we 

use 

equat 1 on 

13 to 

obtain 

the 


dependence of on the local abundance ratios 

and pilfi. . In so doing we have assumed an energy of 
100 MeV/nuc, have included the effects of ionization 
energy loss and have used nominal values of cross 
sections, spectral indicies and elemental abundances of 
parent species. However, we have individually randomized 
the local isotopic composition of each parent element in 
order to avoid biasing the results toward a solar-like 
abundance distribution. As shown ip section 3 3*3*i 
systematic shifts of the calculated source abundances can 
be produced by varying one's assumption concerning the 
local isotopic composition of the parent elements For 
this reason we note that if observational information 
concerning the isotopic make-up of the dominant parent 
elements Is available it may be inappropriate to use the 
curves which we will present when attempting to derive 
source ratios from cosmic ray abundance measurements 
Instead, these curves should be recalculated incorporating 
as many data as are available regarding the isotopic 
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composition of the parents. However, the curves which we 
present are useful for investigating the dependence of the 
various uncertainties discussed in section 3>3-2 of +he 
source and local ratios. 

The other two families of curves which we present 
show the dependence on the local abundance ratios of the 
dominant uncertainties in the derived source ratios. 
These uncertainties are due to statistical errors in, the 
measurement of the local abundances and to uncertainties 
in the parameters in equation The calculation of the 
statistical uncertainties extends the results presented in 
section 3-3»2*5 for the case of a solar-like cosmic ray 
source. In evaluating the calculation uncertainties we 
restrict our attention to the effects of uncorrelated 
errors in the fragmentation cross sections. At the 
present level of uncertainties the effect of these errors 
significantly exceeds the combined effects of the other 
errors discussed above. However these other errors set a 
lower limit on the level to which the uncorrelated cross 
section errors can be reduced before the curves which we 
present can no longer be interpreted as representing the 
overall ca I cu 1 at i ona I uncertainty in obtaining q./q,. 

•I ** 
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3*3-3*i Source Ratio Curves 


The curves relating the source ratio q*/qi, to the 
local abundance ratios and are obtained in the 

following manner. Using nominal values of cross sections, 
spectral indicies and elemental abundances of parent 
species, but randomly varying the isotopic composition of 
each parent element, we have evaluated equation a 
number of times for each of a large number of combinations 
of the observed ratios, PJpy. and . For each such 

pair of local isotope ratios, the resulting values of the 
source ratio, qj/qj^^» were averaged to obtain an unweighted 
mean of those source ratios which are consistent with the 
observed elemental abundances of heavier species. For 
selected values of the source ratio, q^/qj^, interpolation 
has been performed to obtain the locus of points in the 
vs. qj/qj< plane which can be obtained from a cosmic 
ray source in which the species i is absent and which 
yields the desired local ratio, 9hlp}f, • after propagation. 
In addition, these combinations of source and local 
abundance ratios are consistent with a leaky box 
propagation model and with observed values of local 
spectra and elemental abundances of parent nuclei. These 
curves are plotted as solid lines in figure 32a through e. 
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FIGURE 32 


Families of curves showing the relationship 
between cosmic ray abundance ratios in local 
interstellar space and abundance ratios at the 
sources (solid curves). Also shown are the 
uncertainties which result when source ratios 
are calculated from observed abundances due to 
1) statistical errors in the measured local 
abundances (dotted curves) and 2) errors in 
calculating the secondary corrections to the 
observed fluxes due to uncertainties in the 
nuclear fragmentation cross sections (dashed 
curves). Calculations are based on the 
secondary tracer technique (see text) as applied 
to triplets of isotopes (k, j, i) of a single 
element. Species i is assumed to be absent at 
the cosmic ray sources. The ordinate is the 
source abundance ratio, while the 

abscissa is the corresponding local ratio, 
The solid curves are parameterized by 
the local ratio which characterizes the 
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amount of secondary production during 
propagation. The dotted curves, showing the 
source ratio uncertainties due to statistical 
errors in the flux measurements are based on the 
observation of 1000 events distributed among the 
three isotopes of interest. The indicated 
fractional uncertainties can be scaled in 
proportion to V1000 /n' for observations of N 
events. The dashed curves, showing 
uncertainties due to uncorrelated fragmentation 
cross section errors, are based on cross section 
errors of 35^- The indicated fractional 
uncertainties can be scaled in proportion to 
f/Q*35 for relative cross section errors of 
(lOOxf )?. 

In each case we indicate the point 
corresponding to cosmic rays obtained from a 
source having a solar-1 ike isotopic composition, 
as given by Cameron (1S13). See table 3- In 
cases (c) and (d) we show the abundances which 
would result from production of the elements S 
and Ca, respectively, by explosive oxygen 
burning in environments having various degrees 
of neutron enrichment, 7 ^, as presented by 

(1113). See section 3*3.H of 


Woos ley et a I . 



ns 


the text. 

The isotopes (k, j, i) treated in the 
various plots are as follows: 


a) 


'«0, 

'^0 

b) 

Ne , 


*'Ne 

c ) 



”s 

d) 



‘^^Ca 

e ) 

^^'’ca. 

‘^^Ca, 
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LOCAL '®0/'®0 



SOURCE 22Ne/20Ne 
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LOCAL 22Ne/20Ne 



SOURCE 345/325 
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A number 

of 

features 

are 

common 

to all 

of these 

figures. F i rst. 

not 1 ce 

that 

along 

the line 

we 

have 

11 

qj/q 

. . This 1 
k 

s because 

of our 

assumpt 1 on 

that 

spec 1 es 

1 

IS absent 

at 

the cosmic ray sources. If 

species i is 

a 1 so 

absent i n 

the 

local 

cosmic 

rays, then 

the 

amount 

of 

1 nter ste 1 1 ar 

mater i a 1 

traversed during 

propagation must 

be neg 1 i g i 

b 1 e 

and the 

source 

abundances 


will be the same as the observed abundances. 

Second, at a fixed value of the source ratio, q./q , 
the slopes of the various curves of constant 
increases as one goes to larger values of This 

increasing sensitivity of the local ratio 

results from the fact that as 0. increases relative to 
the fraction of the flux 0^ which is of secondary origin 
also increases. Therefore a given absolute change in 0^ 
will correspond to an increasingly large relative change 

in q. as one increases 

4 

Third, for each value of 0\/0y^ there is a minimum 
value of' which can be attained ( correspond i ng to 

qVq^ = 0>- If the measured ratio 0^/0,^ falls far enough 
below this limiting value of the curve correspond i ng to 
the measured ratio 0i/0f^ to be inconsistent with errors in 
measuring these ratios and in calculating the location of 
the curves, we would conclude that species i is unsuitable 
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for use as a tracer of secondary production and probably 
has a non-neg I i g I b I e source abundance. Note, however, 
that it is possible for a point to fall above the limiting 
value of even if q.>Q. in some cases, for example 
the isotopes of Ar, it should be possible to use 
consistency checks between source abundances obtained 
using more than one possible tracer to eliminate this 
poss I b 1 I 1 ty . 

8*3-3*8 Statistical Error Curves 

In calculating statistical uncertainties we have 
assumed that we have 1000 events distributed among the 
three isotopes of interest. The numbers of events of each 
of these isotopes are assumed to be statistically 
independent. The uncertainty in the deduced ratio can be 
expressed as 
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The partial derivatives of with respect to the 
measured ratios and are approximated by finite 
differences obtained from the curves in figure 38- 

The dotted lines in figures 32a through 32® a>“® lines 
of constant fractional uncertainty in the calculated 
values of statistical errors in the 
measurement of local abundances. These curves are labeled 
with the fractional error obtained when a total of iOOQ 
events are accumulated. The uncertainties can be scaled 
by a factor of '\fl000/i^ if +he total number of events is 
N. 

3.3*3*3 Ca I cu I at I ona I Error Curves 

In order to evaluate the uncertainty in source ratios 
caused by the uncertainty in the location of the various 
curves of constant qj/q. m figures 32® through 32® we 
have repeated the procedure used to obtain these curves 
while randomly varying the fragmentation cross sections 
and the isotopic abundances of the parent nuclei. As 
shown in section 3-3*2 the uncertainty in the calculated 
source abundances is due almost entirely to uncorrelated 
errors in the fragmentation cross sections, if the local 
abundances are precisely known. Therefore, this 
calculation of the uncertainties should closely 
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approximate the overall uncertainties which would be 

obtained due to the combined effect of the errors in al I 

the parameters discussed in section 3-3-8> other than the 

flux ratios and calculating the rms 

deviation of a number of values of generated at 

constant 3nd we obtain the required 

uncertainties. In figure 32 we indicate b.y dashed lines 

the curves along which various fractional uncertainties in 

q./q. result due to the uncertainties in calculating the 
i K 

locations of the lines of constant q./q . The qualitative 

J ” 

characteristics of these uncertainty lines are easily 
understood. The uncertainty in q./q, becomes large at 

J ^ 

large values of independent of the value of 

since when becomes comparable to ^ the flux ^ must be 
largely of secondary origin and hence, the rather large 
uncertainties in secondary production will be directly 
reflected in the deduced source ratio. When P; is small 
compared to p. we only get a large uncertainty in qj/q. 
when IS also small. Again this simply means that 0^ 

has a large secondary component (although does not) and 
the uncertainty in this component will be reflected in the 
uncertainty in the deduced source abundance. 


The lines of constant percentage error due to 
uncorrelated uncertainties in the fragmentation cross 
sections are based on an rms uncertainty of 35^ these 



cross sections. 


The calculation uncertainties rndicated 


in figure 32 can be scaled in proportion to this 
percentage uncertainty. 

In the case of Ne (figure 32^^ 'fee local ^*Ne/^Ne 
ratios in excess of '^Q.35> +he ems spread (due to 35^ 
cross section uncertainties) in the source value of 
■*’’^e/^‘’Ne becomes larger than would be obtained by scaling 
from the case of cross section uncertainties smaller than 
35/S* This IS because in this region the secondary 
component of the local flux of ^‘^Ne Is becoming a 
significant part of the total flux, so the denominator in 
equation 12 is relatively small, The 35/^ variation of the 
fragmentation cross sections produces a distribution of 
source ^e/^®Ne values having a long tail. This tail 
causes an unusual ly large value of the rms spread of these 
source ratios. Since we are attempting to produce a 
family of curves which can be scaled to yield source ratio 
uncertainties due to cross section errors we have 
recalculated the curves In the region where the local 
value of ^'Ne/*‘'’Ne is greater than Q.35 using cross 
section errors of il.5?* These results (after being 
scaled up. by a factor of 2^ plotted in figure 32^ 
rather than the corres pon d i ng results obtained for cross 
section errors of 35^* In this way we obtain a family of 
curves which can be scaled to smaller cross section errors 
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over the entire range of local abundances plotted in the 
f 1 gure . 

A similar situation occurs in the case of S (figure 
32c) for local ratios exceeding The 
curves have been modified in this region in the same way 
as described above for the case of Ne None of the other 
three cases considered required such a correction over the 
range of parameters being considered. 

Note, however, that when uncorrelated cross section 
errors are reduced to ;$10^ the other contributions to the 
uncertainty in calculating the location of these curves 
become significant and these curves can no longer be used 
to obtain the overall uncertainty in the calculation of 
source abundance ratios. 


3-3-H Range of Source Abundances 


In the preceding sections we have shown that over a 
sizeable range of local abundances of the isotopes of []• 
Ne, S and Ca large uncertainties will result when source 
abundance ratios of these isotopes are calculated both duo 
to the poor statistical accuracy of the measurements which 
are presently possible and to the large uncerta i nt i es in 
cross sections required for estimating the secondary 





component in the observed fluxes. It has also been shown 
(see table 8) that for soiar-ltke source abundance 
distributions of the isotopes of these elements such large 
uncertainties will be obtained in most of the cases 
considered here. We now wish to examine the question of 
whether the errors expected if the cosmic ray source 
composition is of some reasonable, non-solar form may be 
less significant. Such an examination will aliow us to 
determine the levels to which the errors considered in 
section 3*3*3 must be reduced in order to yield 
significant results. 


We have considered the range of abundances obtainable 
in the charge range if cosmic rays in this 
charge interval are produced by explosive oxygen burning 
(Woosley et al., iS13)* These authors show that for 
appropriate choices of the temperature, density and 
neutron excess In the explosion it is possible to 
accurately reproduce the observed solar system abundances 
in this charge range. They also examine the effect of 
varying the neutron excess on the resulting abundances. 
The neutron excess is defined as 



'n„-ni. 

■nv>+T)p 


where and n^ are respectively the number densities of 



neutrons and protons (both free and bound tn nuclei) in 
the region of the explosion. Solar-1 ike abundances are 
produced for • For selected values of T^, by 
using the calculated source ratios **^Ca/*^^Ca, ^^Ar/^^Ar, 
^’^CI/^*CI and given by Woosley et al. <1S13) and 
assuming that the other isotopes of these elements are 
absent at the cosmic ray sources we have calculated the 
expected local abundances resulting from propagation with 
an exponential distribution of path-lengths with a mean of 
5.5 g/cm*. We have performed the calculations at 
100 MeV/nuc and have neglected ionization energy loss. 
The source abundances used for isotopes of the elements 
other than Ca, Ar, Cl and S were those of Cameron (1013) 
as described in section 3*3-3- 

-2, - 1 

In varying the neutron excess from 10 to If) a 
hundred fold increase in the source ratios S/ S, 
■^*^Ar/^*Ar and *^®Ca is produced. In figures 32^ and d 
we show the points at which the local abundances resulting 
from this calculation fall for selected values of 7^, The 
results for the Ar isotopes are not shown due to 
difficulties in interpreting Ar abundances without 
information on Ca abundances, as discussed above. It can 
be seen in these figures that an increase of by a 
factor 3 over the value which produces solar abundances 
would yield source abundances which should be marginally 
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FIGU RE 22 


Numbers of events (solid curves) and level of 
fragmentation cross section errors (dashed 
curves) required in order to distinguish cosmic 
ray abundances synthesized by explosive oxygen 
burning in environments with various levels of 
the neutron excess, , from abundances 
synthesized in an environment whose composition 
is solar-like. Also shown is a scale indicating 
the collection factors required in order to 
obtain various numbers of events in observations 
made near 500 MeV/nuc. The values of at 
which the explosive oxygen burning model 
reproduces the solar value of the Isotope ratios 
being considered are indicated. Finally, in the 
case of the calcium isotopes, we indicate by an 
arrow the level of fragmentation cross section 
errors below which other sources of 
ca I cu I at I on a I uncertainty exceed the uncertainty 
arising from the cross section errors (in the 
case of sulfur, the cross section errors are 
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dominant over the range 

which 

we are 

cons I der 1 ng ) . 

The triplets of 

1 sotopes 

emp 1 oyed 

in obtaining 

the curves are: 



a) 




b ) a , 

■•"ca, ’^Ca. 





NUMBER OF EVENTS REQUIRED 


Rl 



original 
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FRAGMENTATION CROSS SECTION ERRORS 









distinguishable form solar-like abundances in the case of 
S and indistinguishable fn the case of Ca 


We have used the dependence of the isotopic ratios 

and '^^Ca/ '^‘’Ca on from Woosley et al. (1^113) 

the solar values of these ratios from Cameron (1913) a>^d 

the error calculations described above to estimate the 

degree to which measurements of local fluxes and 

measurements of fragmentation cross sections must be 

improved in order to distinguish between production in a 

solar-like environment and production in a region with 

enhanced neutron excess. We have assumed that the two 

sources of error (measurement statistics and cross section 

errors) contribute equal ly to the error in determining 

q /q and we require that the overal I error in qj/q. be 
J K J * 

50^ of the difference between the values of 

solar case and for the value of yj being considered. 


In figure 33a we show, as a function of j the 
number of sulfur events (solid curve) needed in order to 
reduce the statistical uncertainty in the calculated value 
of the source ratio far enough to distinguish this 
ratio from the correspond I ng solar-like ratio. The 
collection factor ( m^sr sec ((^eV/nuc)) required in order 
to obtain these numbers of events above the earth's 
atmosphere at an energy /^5H0 MeV/nuc is also shown The 
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dashed line shows the level to which uncorrelated errors 
in measurements of the relevant fragmentation cross 
sections must be reduced in order to make the same 
distinction between solar and non-solar environments. 
Figure 33b shows the correspond i hg values for the case of 
the source ratio a / a . In this case the calculated 
fractional error in the fragmentation cross sections must 
be reduced to a level which is so low that other 
calculation errors will no longer be negligible. In table 
8 it can be seen that when the uncorrelated fragmentation 
cross section errors are reduced to <^1^ (in the 
solar-like case) they will contribute an amount equal to 
the combined effect of al I other calculation errors to the 


error in q^/qj^. Therefore even if the fragmentation cross 
section errors can be reduced to ^10?, it will not be 
possible to distinguish a cosmic ray source with ?^^<4xl0 
from a solar-1 ike source unless other errors discussed in 


section 3*3-2 can also be reduced. In figures 33a and b 


we indicate the value of at which the model of Woos I ey 
et al. (1913^ reproduces the solar system value of the 


abundance ratio being considered. The difference of the 


two values of yj is well within the accuracy expected from 
the nucleosynthesis calculations. 





3 . 3*5 Solar Modulation Effects 

In the above discussion we have assumed that 
observations could be made of abundances of the cosmic 
rays in local interstellar space. However, fluxes 
observed near earth have been modified by their 
interaction with the fluctuating magnetic field in the 
solar wind in penetrating to the orbit of earth. This 
phenomenon of solar modulation has been extensively 
studied. The subject has been reviewed by Jokipii 
and by Fisk (1514). 

The effects of convection, diffusion and adiabatic 
energy loss on cosmic ray spectra are described by a 
Fokker-P i anck equation (Parker, 14b5). These effects will 
alter spectra by preventing a fraction (depending on 
magnetic rigidity) of the interstellar particles from 
penetrating to the earth's orbit and by redistributing 
particles having a particular energy in interstellar space 
over a range of degraded energies. 

Gleeson and co-workers have shown (Urch and Gleeson, 
1513; Webb and Gleeson, 1513; Gleeson and Webb, 1515) 
that those particles with a mass-to-charge ratio of 3 
which are observed at a particular energy near earth 
originate in interstellar space from a distribution of 
energies whose mean is 'wl0(]-150 MeV/nuc greater than the 



observed energy during times of minimum solar activity and 
‘“200~350 MeV/nuc greater than the observed energy at solar 
maximum. These distributions are approximately Gaussian 
in shape with ful I width at half maximum in the range 
30~50 MeV/nuc at solar minimum and ~|0-ii0 MeV/nuc at solar 
maximum. In addition, particles of lower M/Z but the same 
energy per nucleon lose more energy during solar 
modulation than those of higher M/Z The difference in 
the mean energy loss by particles with M/Z=2*E and those 
with M/Z=2 IS *^20 MeV/nuc at solar minimum and MeV/nuc 
at solar max imum . 

Isotopic abundance ratios measured at earth will 
differ from those which are present at the same energy in 
interstel I ar space due to two effects. First, due to the 
loss of energy during solar modulation the particles at a 
fixed energy at earth originated at higher energies in 
interstellar space. If the ratio being considered is 
strongly dependent on energy then the magnitude of the 
energy loss must be accurately determined in order to 
interpret the measured abundances. This situation occurs 
when one is considering the ratio of an isotope which has 
a large secondary component to one which is largely 
primary if observations are made below a few hundred 
MeV/nuc, since at these low energies many of the 
fragmentation cross sections exhibit a sizeable energy 
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dependence and the effects of ionization energy loss on 
the spectral shape become significant. At energies 
approaching iOOO MeV/nuc, however, the fragmentation cross 
sec-tions are approximately independent of energy and 
ionization energy loss effects are unimportant so 
abundance ratios should approach a constant value. 

Second, even if the interstellar spectra of the 
species of interest have a constant ratio at al I energies, 
the ratio measured at a constant energy per nucleon at 
earth will be altered since different isotopes with the 
same energy per nucleon have different magnetic rigidities 
(proportional to their mass-to-charge ratios). The 
transmission of interstellar particles to earth and the 
redistribution of their energies is predom i nant I y a 
rigidity-dependent process since it arises from the 
scattering of the particles from Irregularities in the 
interplanetary magnetic field. 

In order to estimate the magnitude of this effect we 
have performed a numerical solution of the spherically 
symmetric Fokker-P I anck equation 

0=ji5[r“K|“ (15) 

using the method of Fisk In this equation r is 
the radial distance from the sun, T is the kinetic energy 
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IS the solar wind 


of the particles being considered, 
velocity (taken to be 400 km/sec), K Is the diffusion 
coefficient. In addition, U is the number of cosmic ray 


parti c 1 es 

per 

unit 

VO I ume and 

per 

unit kinetic 

energy 

( re 1 ated 

to 

the 

flux, j , 

by 

j=^cU/4‘ii') 


and 

OC = ( T+2mc 

^)/(Ttmc^) 

(mc^ being 

the 

rest energy 

Of 

the 

part Ides 

being 

cons 1 dered ) . We 

have 

assumed that 

the 


interstellar spectra of both species of interest are 
proportional to the interstellar spectra obtained by 
Garrard (1413^ from studies of the modulation of H and He 
nuclei : 

j (m 0.75-(?3l Me\^/nuc)) 


The diffusion coefficient has been assumed to 
form 


K 



Pa/? 

P<P. 


be 


of 


the 


where 

and 

are 

constants. 

P IS the 

part i c 1 e 

»s 

magnet i c 

rigidity 

and 

^c is its 

velocity. 

Th 1 s form 

1 s 

assumed to 

apply w 1 

1 th i n 

a radius D 

of the sun 

and it 

1 s 


assumed that there is no modulation outside of this 
radius. We have used Garrard's values for P and for the 
modulation parameter 
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7^(Mu)= (,5) 

where is the solar wind velocity and is assumed to 

have a constant value of HOO km/sec. Using values of 
D = J.O, 20 snd 50 AU we have used equation |5 to calculate 
k^. The parameters k<j, and D along with the spectrum 
15 were used in obtaining spectra at 1 AU from equation 
13. In table 5 we show the abundance ratios obtained at 
earth at a fixed energy per nucleon as a function of this 
energy for species having mass~to-charge ratios of 2 
2.2- lb addition to the results obtained from the 
numerical solution using the parameters given above, we 
show results obtained using the force field approximation 
(Gleeson and Axford, i%8) with the same values of 7|( lAU) 
given by Garrard. It can be seen from the table that the 
effects of solar modulation are not sensitive to D if the 
value of yjC lAU) ts held constant. At solar minimum 
direct comparison of abundances measured at a fixed energy 
at 1 AU will result in an error of rv5-15^ when comparing 
isotopes with M/Z=2 sbd M/Z = 2-2 ^^ue to differences in the 
level of modulation at the different rigidities of these 
particles. At solar maximum this error increases to 
'^iO'25/S. It will be necessary to apply the appropriate 
solar modulation corrections to abundance rat los jmeasured 
near earth if solar modulation effects are not to be a 



Tab 1 e £ 

Correction Factors for Solar Modulation 


E/M 

— 0 = 

10 AU — 

— D = 

20 AU — 

-D=50 

AU -- 

force 

field 

( MeV/nuc ) 

i0b5 

1410 

14b5 

1410 

IM 

14_1Q 

l_4b5_ 

mo 

100 

M3b 

1.233 

M33 

1-228 

M43 

1-225 

1.041 

1-208 

205 

MOM 

1.20b 

M02 

1.203 

M30 

1.203 

1.011 

M81 

31b 

1.088 

1.118 

1.084 

Mil 

1.088 

M14 

1.0b4 

Mb3 

5b2 

l.Oba 

M40 

1.0b2 

M40 

1.054 

1 153 

l-05b 

M3b 

150 

1.041 

M21 

1.052 

M23 

\ 

1.05b 

M43 

1-050 

M21 

1000 

1.045 

1.102 

1.041 

MOb 

1-045 

MU 

1.043 

MOb 


* Cj(M/Z=2.^)/j(M/Z=?)]^^^ /Cj(M/Z=B.?)/j(M/Z=p)D^5 
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serious limi+atlon on our ability to interpret such 
measurements . 


Finally, we note that with the limited counting 
statistics achievable with present day instruments ft will 
be necessary to combine data over a fairly large energy 
interval (at least several hundred MeV/nuc) in order to 
reduce statistical errors to an acceptable level. As a 
consequence, the averaging resulting from the solar 
modulation process should not presently severely limit the 
capability of these experiments. 
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4. SUMMARY AND CONCLUSIONS 


We have presented a detailed evaluation of the mass 
resolution capabilities of the Caltech HEIST instrument. 
The instrument employs an energy~!oss - residual energy 
technique to make mass measurements of elements from Li 
through Ni at energies between 30 and 8QQ MeV/nuc. We 
have shown that Landau fluctuations in the energy loss 
measurement result in a fundamental limit to the 

achievable rms mass resolution which ranges from '^♦01 AMU 
to <^.2 AMU (Z'v^b). Additional contributions to 
the mass uncertainty due to errors in measuring the 
trajectory and energy losses of cosmic ray particles have 
been shown to increase the overal I rms mass resolution to 
i AMU for and to 3 AMU for Z~2b- Such 

resolution will permit unambiguous separation of adjacent 
isotopes of the lighter elements and separation of 

isotopes differing by two mass units in the iron region 
(assuming that the isotope lying between these two has 

negligible abundance). 


Cosmic ray isotope spectrometers will, within a few 
years, have provided data which will make it possible to 
tell whether the isotopic composition of the cosmic ray 
sources is dramatically different from the composition of 
solar system material. For example, if the composition of 
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an element in the cosmic ray sources is dominated by a 
single isotope different form the isotope dominating the 
solar abundances, this difference should be readily 
identified. 

We have gone on to evaluate in detail the 
uncertainties which will arise in deriving abundance 
ratios at the cosmic ray source from measurements made in 
local interstellar space. In this analysis we have 
considered elements which have at least one isotope which 
IS likely to have a negligible abundance at the cosmic ray 
sources. This element is employed as a tracer of the 
production of secondary cosmic rays by means of spallation 
reactions during propagation. By basing secondary 
corrections to the observed fluxes on the observed 
abundance of a purely secondary isotope of the same 
element, we are able to avoid uncertainties resulting from 
the possibility of charge dependent acceleration 
mechanisms. This analysis was performed in the context of 
the leaky box model of cosmic ray propagation. 

The uncertainties affecting the calculation of cosmic 
ray source ratios are of two types: 1) those arising from 
the statistical errors in the measured fluxes and 2) those 
arising from uncertainties In the calculation of the 
secondary component of these fluxes. In typical present 



day cosmic ray experiments the statistical uncertainties 
in the der i ved .sour ce abundances should be significant, in 
some cases exceeding the abundances themselves. 

The uncertainty in calculating the secondary 
corrections arises from the combined effects of 
uncertainties in a number of parameters employed in this 
calculation. The dominant contribution arises from 
uncorrelated errors in the important nuclear fragmentation 
cross sections. Many of the important cross sections have 
not yet been measured and are typically obtained from 
sem I -emp I r I ca I formulas. The uncertainties in these cross 
sect { ons~~est I mated to be r.' 35/5~~a I so produce source 

abundance uncertainties comparable to the abundances 
themselves. Measurements of the important cross sections 
can, of course, reduce this uncertainty. We have shown 
that it should be particularly useful to measure with high 
precision the ratios of cross sections for producing the 
isotopes of interest from individual parents. Even if 
errors ^^35^ remain in the absolute cross sections, the 
effect of these errors on the calculated source ratio will 
be reduced by a factor ^3- 

Other calculation uncertainties will become 

N 

significant when the fragmentation cross section errors 
can be reduced by a factor of 3 or 4 from the 35^ level 
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which we have assumed. The most important of these are 
due to: |) uncertainties in the composition (both 

elemental and iso+opic) of the parent species which are 
fragmented to produce the observed secondaries and 
uncertainties in the shape of the equilibrium interstellar 
spectra of the isotopes which we are considering. 

By assuming that all stable isotopes of a given 
element are, a priori, equally likely in the local cosmic 
rays, we obtain uncertainties of ;$ 50 ^ 'h the calculated 
source ratios. In addition, we have shown that it is 

possible with specific parent isotope distributions to 
obtain significant systematic shifts of the calculated 
source ratio from the value obtained by averaging over 
randomly assigned isotopic compositions. For cosmic rays 
derived from a sotar-l ike source, these shifts can be 
several , times the standard deviation obtained from 

randomizing isotopic compositions. We have pointed out 
that since often the production of secondary cosmic rays 
is dominated by the spallation of species of approximately 
the same atomic number, it is likely that simultaneous 
measurements of the isotopic composition of an element of 
interest and of a few important parent elements will 
substantially reduce both the rms error and the systematic 
shifts which we obtain from randomized isotopic make-up of 
the parent element. 
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Uncertainties in the elemental composition of the 
parents also yield a source ratio uncertainty of ;S50?- 
Our estimates of these uncertainties are based on 
elemental composition errors which largely reflect the 
variation from experiment to experiment of the measured 
elemental abundances. These d i f f erences--presumab ! y due 
to systematic errors in the measurements — should soon be 
reduced by measurements using higher resolution 
instruments with' good statistical accuracy. Therefore the 
present uncerta i nt i es in the elemental composition of the 
parent population should not seriously limit our ability 
to calculate accurate source abundances. 


Uncertainties in the shape of the spectra of the 
species of interest contribute source ratio errors which 
are typically ;S30^* Precise determination of the 
interstellar spectra requires flux measurements with large 
statistics at a number of energies and also accurate 
calculations of the effects of solar modulation. It is 
clear that such refinements will not be available for 
interpreting the data from the first generation of cosmic 
ray Isotope experiments. We note, however, that the 
effects of uncertainties in the spectra will be minimized 
for measurements at energies in excess of [ GeV/nuc since 
at these high energies the effects of ionization energy 
loss are minimized (and it is these effects which make the 
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spectral shape important). Also, at high energies the 
influence of solar modulation on the observed spectra will 
be less than at lower energies. 

Families of curves were presented for the elements 0 , 
Ne, S and Ca (figure 32) which show the dependence of the 
resulting source ratio uncertainties on the observed 
fluxes. In computing the effect of the uncertainty in the 
calculation of the secondary contributions we have" only 

i 

considered the errors in the nuclear fragmentation cross 
sections, since presently these constitute by far the 
dominant source of ca I cu I at i ona I error. Curves showing 
the statistical errors were based on a total of 1000 
events distributed among the three isotopes of the element 


being considered. The 

curves showing the 

errors due 

to 

fragmentation cross 

sect 1 on 

errors 

were based 

on 

uncorrelated errors of 
! 

1 h these cross 

sect ions. 

The 

scaling of both types 

of errors 

has been d 

1 scussed . 


Using these curves we have 

exam i ned 

the extent 

to 


which the uncertainties must be reduced in order to 
distinguish fairly minor differences In the isotopic 
composition at the cosmic ray sources. In so doing we 
have special ized to the case of abundances produced by 
means of explosive oxygen burning in environments 
exhibiting various levels of neutron excess. Using the 
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neutron excess as an independent variable we have shown 
the number of events and the level of fragmentation cross 
section errors which will be required in order to 
distinguish production in an environment with this neutron 
excess from production in a solar-like environment. 

Finally we have shown that solar modulation effects 
are capable of altering the abundance ratios observed near 
earth by up to 25^ from those present outside the 
heliosphere. Consequently solar modulation corrections 
must be made as a first step in interpreting any isotopic 
abundance observations made near earth. 
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APPENDIX A -- UNCERTAINTY IN MEASURED ABUNDANCE RATIOS 

In this appendix we obtain estimates of the 
uncertainty in abundance ratios obtained from measured 
mass distributions. In section A.| we obtai'n a 
relationship between the desired uncertainty in a ratio 
and the number of events required in order to reduce the 
uncertainty to this level. This relationship depends on 
the magnitude of the ratio and on the mass resolution. In 
section A. 2 we summarize the derivation of the relation 
between mass resolution and abundance ratio required to 
obtain an inflection point in the measured mass 
distribution. This relationship is a measure of the mass 
resolution required in order to insure that systematic 
errors cannot invalidate the deduced abundance ratios. In 
both of these treatments we consider only the case of two 
isotopes with non-neg 1 I g i b I e abundance and we assume that 
the mass resolution is equal for the two isotopes. 

A.i Statistical Errors 

We consider a mass distribution consisting of two 
Gaussian mass peaks having unit separation between their 
means, standard deviation CP and relative abundance r. 
The probability density for samples obtained from this 


d I str I b ut i on is: 
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where f=r/(l+r). If a sample of N events is obtained from 
this distribution, one can obtain an estimate of f by 
various means. One technique is to accumulate the 
measured events in mass bins and to use the maximum 
likelihood method (see, for example, Mathews and Walker, 
iSlQ) to obtain the most probable estimate of f. We wilj 
represent the number of counts in the i-~ bin by n. and 
the lower and upper limits of the i bin by M, and M,^, 
respectively. Then the probability of any particular 
event falling in the bin can be written as 



oo 


Also, we denote , by N the total number of events 
accumulated and by ^ the uncertainty in the number of 
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events obtained in the bin. 

[f the number of counts in each bin has a Gaussian 
distribution, then several simplifications result. First, 
the maximum likelihood estimate of f is simply the value 
of f obtained by performing a least squares fit of A.l to 
the data: 



Second, the uncertainty in our estimate of f can be 
obtained in simple analytical form as 


I 



In the case which we are considering, the number of counts 
in each bin will have a Poisson distribution with mean Np; 
and, hence, standard deviation However the central 

limit theorem of statistics (Mathews and Walker, 1410) 
guarantees that for sufficiently large values of Np* this 
distribution approaches a Gaussian with these values of 
the mean and standard deviation. For the sake of 
simpi I c I ty we employ equation A. 3 in al I cases using 
to obta I n 


I 
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The simple scaling of 0^ in proportion to I/'ifw’ is likely 
to breakdown for smalt values of N. 

The uncertainty in the abundance ratio, r, can be 
obtained from ^ using 

- OF 

We have evaluated equations A.H and A. 5 using 3b bins 
each of width AM=0*2 between =H order 

to obtain the curves shown in figure 3* have, found 

that the values of 0^ obtained vary by less than 10^ when 
AM IS varied between 0-Q5 and 0*5 and when is varied 

between -3 and -1 (with +i) • 




A. 3 Systematic Errors 

As discussed in section 3*1 "bbe possible presence of 
systematic errors in the determination of particle masses 
makes it important to achieve sufficient mass resolution 
to obtain separate mass peaks for adjacent isotopes. As a 
criterion for separation we demand that Gaussian 
distributions with means separated by one unit have an 
inflection point. Mass resolution which is better than 
the value obtained by applying this condition will yield 
two distinct maxima in the mass distribution whereas worse 
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resolution will not. In the following treatment we will 
assume that no statistical uncertainties are present 
(i.e., that a large number of events has been 
accumu I ated ) . 

We are considering the distribution p(M) given in 
equation A.l. The conditions for obtaining an inflection 
point in this distribution are p*(M) = Q and p''(M)=[|. 
Applying these conditions to equation A.i we find that for 
f<i/2 there will be an inflection point at 

M=i 

if O" and r=f/(l-f) are related by 



This relationship between the mass resolution and the 
abundance ratio is shown in figure H- 
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APPENDIX B — TRAJECTORY UNCERTAINTIES . 

In this appendix we evaluate the uncertainties in the 

trajectory of a particle as determined by measurements in 

our proportional counter hodoscope. In particular, we 

obtain uncertainties in the position at which the 

trajectory intersects various levels in the scintillator 

1 

stack and the uncertainty in sec6, where 6 is the 

\ 

particle's angle of incidence measured from the axis of 
the crystal stack. 

We consider the case of four measurements of the 
x-coordinate and four of the y-coordinate made at various 
positions along the particle's track, and we assume that 
each of these measurements has an uncertainty of 
We illustrate this geometry in figure B.i. The origin of 
the z axis is chosen to coincide with the depth at which 
we wish to obtain the uncertainties in the absolute x and 
y positions of the track. We separately fit straight 
lines to the x and to the y measurements: 
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F IGURE B.l 


Schematic illustration of parameters 
used to describe the hodoscope geometry 
(Zo» trajectory of a 
particle (x^, dx/dz). Values of Az, and 
Az^ for the HEIST hodoscope are shown. 



Az| = 4 48cm 
Az 2 = 30.61cm 



XI MWPC 
X2MWPC 


X3MWPC 
X4 MWPC 


ru 

h— 


X 



2?0 


The maxttnum likelihood method {Mathews and Walker, iSIO) 
can be used to evaluate the parameters x^, , dx/dz, and 
dy/dz. In this case the maximum likelihood method reduces 
to a simple linear least squares fit and the coefficients 
are obtained by solving the matrix equation 



The uncertainties in the fitted coefficients are 
obtained from the diagonal elements of the ’’error matrix". 
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^(Az,y-+(Az^^' 


The analogous expressions are obtained for and dy/dz. 
Note, however, that in order to obtain x^ and y at the 
same depth In the crystal stack, different values of z, 
must be used since each y measurement is made at a depth 
which is approximately 2-24 cm closer to the crystal stack 
than the corresponding x measurement. 


B.l Absolute Position Uncertainty 


In table B.l we list the ratio of the position 
uncertainties, and at various depths in the HEIST 

scintillator stack to the position uncertainty, OfTw?c > in 
each of the proportional counters. Values of Az, and 
listed in figure B.l are employed. It can be seen from 
table B.| that this ratio varies between, 0«81 and 1-04 
depending on the depth in the scintillator stack. In our 
geometry the presence of redundant position measurements 
compensates for the degradation of position resolution due 





Position Uncertainties in the Scintillator Stack 


Location in 

Sc 1 n 1 1 1 1 ator Stack 

Zo ■fo'* 

X Measurement 

for 

y Measurement 

0^ 

gTo, 
0 81 

Top of 

DO 

b.1b 

cm 

4.52 

cm 

0.43 

Top of 

D1 

1.10 

cm 

4 8b 

cm 

0.44 

0.88 

Top of 

DO 

1-45 

cm 

5.21 

cm 

0.45 

0.84 

Top of 

D3 

1.11 

cm 

5.53 

cm 

0.4b 

0 40 

Top of 

D4 

8.20 

cm 

5.4b 

cm 

0 41 

0.41 

Top of 

D5 

8.81 

cm 

b.51 

cm 

0.44 

0.43 

Top of 

Db 

4.bl 

cm 

1.31 

cm 

1.01 

0.45 

Top of 

D1 

10.82 

cm 

8.58 

cm 

1.04 

0.48 

Bottom 

of Dl 

12.52 

cm 

10.28 

cm 

1 04 

1.03 
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to the extrapo I at I on of the trajectory from the hodoscope 
to the scinti I lator stack. 


Uncertainty tn sec 6 


We can 
express ions 
that 


obtain the uncertainty 
obtained above for eg* 


I n sec B 

and 0^ 
«<z 


from the 
by noting 


sec 0 



1 


so that 



After simple algebraic manipulations we obtain 


Seed ^ 

Insertion of values of 4z, and from figure B.X in this 

expression yields 


secO 


O.OOIG Sln,3.0-C7~. 


where o;;wfc is expressed 


in millimeters. 
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